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ABSTRACT 


^Design  methods  are  developed  for  determining  aerodynamic  control  effectiveness 
at  hypersonic  speeds  along  body  surface.  Pressure  and  heat  transfer  distri¬ 
butions  in  separated  regions  due  to  aerodynamic  control  deflection  are  des¬ 
cribed  in  terms  of  characteristic  magnitudes  and  distance  parameters  by  semi- 
empirical  correlations.  The  forms  of  these  correlations  are  inferred  from 
theory  and  experimental  data.  Using  these  correlations,  pressure  distribution 
in  the  separated  region  over  a  deflected  flap  is  approximated  and  expressions 
for  force  and  moment  coefficients  are  determined.  General  charts  are  develop¬ 
ed  which  present  separation  and  flow  parameters  over  a  range  of  flight  condi¬ 
tions  for  a  typical  hypersonic  vehicle. 

Flow  separation  over  a  fin-plate  configuration  is  presented  using  experimental 
measurements.  Also,  characteristics  of  flow  over  a  flat  .plate,  flat  delta 
wing,  and  delta  with  dihedral  are  analyzed  using  visual  flow  records  and 
pressure  measurements  from  the  point  of  view  of  two-dimensional  flow.  Appli¬ 
cability  of  the  correlation  expressions  to  separated  f]ow  on  various  config¬ 
urations  is  discussed  and  calculated  aerodynamic  coefficients  are  compared 
with  measured  values.  The  experimental  programs,  test  results,  test  model 
geometry,  correlation  expressions  and  prediction  methods  are  reviewed  crit¬ 
ically,  and  the  applicability  of  prediction  methods  discussed. 
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specific  heat  at  constant  pressure 

lift  coefficient 

pitching  moment  coefficient 

normal  force  coefficient 

constant  in  the  blast-wave  expression 

flap  chord 

drag 

diameter  of  leading  edge 

upstream  interaction  length 

downstream  interaction  length  to  peak  pressure 
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h 

HL 

HM 

K 

K 

k 

L 

L 

L 

4 

4 

M 

M 

N 

n 

n 

n 

Nu 

P 


downstream  interaction  length  to  pressure  rise 
parameter  required  to  describe  boundary-layer  growth 
heat  transfer  coefficient 
hinge  line 
hinge  moment 

hypersonic  similarity  parameter 
general  constant 
thermal  conductivity 
lift 

reference  length  for  moment  coefficients 
length  associated  with  deflected  flap  geometry 
separation  length 
free  interaction  length 

downstream  heat  transfer  interaction  length 

Mach  number 

moment 

normal  force 

coordinate  normal  to  streamline 
coordinate  normal  to  shock  sheet 

parameter  required  to  describe  boundary  layer  growth 

Nusselt  number 

pressure 

total  pressure  behind  the  shock 
pressure  behind  a  reflected  shock 
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8 


Prandtl  number 

dynamic  pressure 

heat  flux 

recovery  factor 

Reynolds  number 

reference  area 

coordinate  along  streamline 

coordinate  along  shock  sheet 

Stanton  number 

temperature 

velocity  of  outer  flow 

velocity  in  viscous  region 

longitudinal  coordinate 

coordinate  of  shock  locus 

distance  to  transition  point 

lateral  coordinate 

vertical  coordinate 

angle  of  attack  or  flow  deflection  angle 
incremental  change 

a  function  of  relative  flap  chord  length 
ratio  of  specific  heats 
flap  deflect j on  angle 
boundary  layer  thickness 

incompressible  boundary  layer  thickness  (no  reference 
temperature  conditions  used) 
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5  wedge  angle 

6*  displacement  thickness  of  boundary  layer 

A  sweep  angle 

H  viscosity 

p  density 

6  shock-wave  angle  measured  from  upstream  flow  direction 

6  inclination  of  surface  to  flow 

<f>  angle  associated  with  geometry  of  separation 

X  hypersonic  interaction  parameter 

(0  exponent  in  viscosity- temperature  relation 


Subscripts  and  Superscripts 


AV 

average 

asymp 

asymptotic 

aw 

adiabatic  wall 

b 

body 

b 

blast  wave 

center  line 

CP 

center  of  pressure 

e 

outer  flow 

fp 

flat  plate 

HL 

hinge  line 

mv 

inviscid 
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incipient 
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SECTION  I 


INTRODUCTION 


The  development  of  efficient  hypersonic  vehicles  involves  solution  of  problems 
associated  with  lower-speed  aircraft  together  with  additional  complexities  in¬ 
troduced  by  the  more  severe  environment  of  hypersonic  flight.  Spatial  arrange¬ 
ment  of  the  vehicle  elements  and  configuration  optimization  for  performance 
and  stability  may  lead  tc  shapes  which  create  complex  flow  fields,  resulting 
in  interference  flow  and  flow  separation. 

Flow  separation  is  a  common  and  important  phenomenon  in  the  aerodynamics  of 
aircraft  and  aerospace  vehicles.  It  can  occur  in  several  ways,  such  as  ahead 
of  deflected  flaps,  on  the  leeward  side  of  a  surface  inclined  at  large  angles 
of  attack,  near  the  impingement  of  a  shock  wave  upon  the  boundary  layer  of  a 
body,  and  on  a  curved  surface.  Separation  often  limits  the  usefulness  of  aero¬ 
dynamic  devices;  for  example,  separation  limits  the  maximum  lift  of  an  airfoil. 

Separation  may  be  referred  to  as  any  reverse  flow  which  increases  drag  or  de¬ 
creases  lift  in  external  flow.  In  low-speed  flight,  the  principal  effect  of 
flow  separation  is  to  cause  drastic  modification  of  the  pressure  distribution. 

In  supersonic  and  hypersonic  flow,  the  presence  of  shock  waves  is  associated 
with  strong  pressure  gradients  which  gives  rise  to  interactions  between  bound¬ 
ary  layers  and  shock  waves.  These  interactions  often  lead  to  flow  separation. 
The  occurrence  of  shock-induced  separation  may  influence  heat  transfer,  pres¬ 
sure  distribution,  and  the  aerodynamic  characteristics  of  a  vehicle.  From  a 
thermodynamic  point  of  view,  the  interaction  region  is  characterized  by  a  re¬ 
duction  of  local  heating  rates  in  the  separated  region  and  a  substantial  heat¬ 
ing  rate  increase  at  the  reattachment  point.  Flow  separation  in  the  vicinty 
of  t railing-edge  flaps  and  around  fins  will  result  in  changes  in  pressure  dis¬ 
tribution  and  will  influence  the  control  characteristics  of  a  vehicle.  It  is, 
therefore,  essential,  from  a  practical  point  of  view,  to  understand  the  flow- 
separation  phenomena  and  to  describe  the  flow  conditions  in  the  separated 
region. 

This  study  is  primarily  concerned  with  the  development  of  semi-empirical  corre¬ 
lations  for  characteristic  parameters  of  a  separated  flow  at  hypersonic  speeds . 
The  correlation  equations  developed  describe  the  pressure  and  heat-transfer 
distributions  in  terms  of  local  flow  properties  by  defining  characteristic  mag¬ 
nitudes  and  distance  parameters.  Knowing  these  quantities,  relations  for  in¬ 
cremental  aerodynamic  force  and  moment  coefficients  are  developed.  Thus, 
methods  are  obtained  for  predicting  the  effectiveness  of  aerodynamic  controls, 
as  influenced  by  separation  phenomena,  over  a  wide  range  of  flight  conditions. 
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Section  II 

FLOW  SEPARATION  FEEHCMENA 


The  fluid  near  the  body  surface  in  a  viscous  flow  is  retarded  by  the  skin 
friction  forces.  The  velocity  of  the  fluid  in  the  boundary  layer  varies 
from  -zero  at  the  vail  to  approximately  the  iaviscid  value  at  the  outer  edge . 
In  the  absence  of  other  forces  acting  to  further  retard  this  flow,  the  slope 
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of  the  velocity  profile, 


is  positive  at  the  vail  and  gradually  ap¬ 
proaches  zero  at  the  edge  of  the  boundary  layer  (Figure  l).  If,  in  addition 
to  vail  shear,  the  flow  encounters  a  pressure  increase  in  the  streamwise  di¬ 
rection,  deceleration  of  the  flow  takes  place-  Due  to  the  vork  done  against 
this  pressure  force,  the  energy  of  the  fluid  is  reduced  and  the  fluid  may  be 
brought  to  rest.  Thj.s  condition  vill  be  experienced  first  by  the  very  low 
momentum  fluid  vithin  the  boundary  layer  and  nearest  the  vail.  Downstream 
of  this  region  the  acting  force  vill  cause  the  fluid  near  the  vail  to  flow 
in  the  upstream  direction,  creating  a  backflow.  The  nonnal  velocity  gradient 

at  the  vail,  —  J  v  must  be  positive  when  the  fluid  next  to  the  vail  moves 
with  the  stream  and  negative  when  the  fluid  in  this  region  flows  against  the 


stream. 


0U  \ 

It  follows  that  where  the  two  flows  meet,  g y/w  = 


This  point 


on  the  vail  which  divides  these  tvo  regions  of  flow  is  defined  as  the  sepa¬ 
ration  point. 


In  most  fluid  dynamic  situations,  the  force  giving  rise  to  separation  is  an 
adverse  pressure  gradient,  i.e.,  the  pressure  increasing  in  the  downstream 
direction.  In  order  for  the  separation  to  occur  the  pressure  rise  must  be 
large  enough  and  the  gradient  must  be  sufficiently  severe.  In  a  supersonic 
flow,  an  adverse  pressure  gradient  may  be  generated  by  the  aerodynamic-  shape 
of  the  body  (for  example,  a  compression  corner)  or  by  an  external  source 
such  as  an  inpinging  shock  wave.  In  each  case,  the  super  sonic  flow  experi¬ 
ences  a  pressure  rise  across  a  shock  wave.  Ideally  the  shock  wave  represents 
a  discontinuous  rise  in  pressure.  Actually  the  pressure  rise  takes  place 
over  a  finite  distance.  There  is  a  region  in  the  part  of  the  boundary  layer 
nearest  the  vail  where  the  flow  is  subsonic.  In  this  region,  pressure  waves 
vill  propagate  upstream  of  the  disturbance.  Experiments  have  shown  that 
such  disturbances  ma>  be  propagated  a  significant  distance  upstream.  The 
process  by  which  this  pressure  field  spreads  through  the  subsonic  portion  of 
the  boundary  layer  is  referred  to  as  "pressure  diffusion." 
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INCIDENT  SHOCK 


An  illustrative  case  of  the  separated  flow  is  an  oblique  shock  incident  on  a 
boundary  layer  along  a  flat  plate. 

An  oblique  shock  impinging  on  a  flat  plate  in  an  inviscid  flow  generates  a 
step  increase  in  pressure  and  causes  a  flow  deflection  as  shown  in  Figure  2. 
When  an  oblique  shock  impinges  on  a  flat  plate  with  a  boundary  layer,  a 
quite  different  flow  pattern  results  (Figure  3)-  The  viscous  flow  cannot 
negotiate  an  abrupt  pressure  rise,  and  therefore  the  pressure  rise  must  be 
spread  over  a  finite  distance.  A  strong  adverse  pressure  gradient  generated 
by  the  incident  shock  causes  the  boundary  layer  to  thicken;  if  the  incident 
shock  is  strong  enough,  the  flow  separates  from  the  wall  and  a  bubble  of 
dead  air  is  created.  Due  to  the  deflection  of  the  external  flow,  caused  by 
the  growth  of  the  boundary  layer  ahead  of  the  shock  impingement,  a  family  of 
compression  wavelets  upstream  of  the  impingement  fans  out.  The  slowly  moving 
air  in  the  separated  region  is  essentially  at  constant  pressure,  and  the  in¬ 
cident  shock  is  reflected  from  it'  in  the  form  of  a  fan  of  expansion  waves. 

The  streamlines  are  deflected  toward  the  plate  and  then  turned  parallel  to 
the  wall,  which  causes  a  flow  compression  and  generates  another  family  of 
compression  wavelets. 

The  pressure  distribution  in  viscous  flow  deviates  markedly  from  that  of  an 
inviscid  flow;  it  shows  a  pressure  rise  ahead  of  the  shock  impingement, 
reaches  a  plateau  value  in  the  separated  region,  rises  to  a  peak  value,  and 
drops  to  the  final  inviscid  pressure  downstream  of  the  reattachment.  If 
the  shock  is  not  strong  enough  to  cause  separation,  the  pressure  pattern  re¬ 
sembles  that  for  a  regular  reflection  except  that  the  rise  is  spread  over  a 
small  length,  of  the  order  of  the  boundary  layer  thickness,  ahead  of  the 
shock  impingement  point.  If  the  shock  strength  is  increased  and  exceeds  that 
needed  for  incipient  separation,  the  separation  point  moves  upstream  of  the 
shock  impingement  point  since  the  boundary  layer  cam  withstand  only  a  weak 
pressure  gradient  in  the  separated  region . 

The  pressure  rise  at  separation  seems  to  be  fixed  by  the  characteristics  of 
the  boundary  layer,  and  to  be  'n~ee  from  any  downstream  influence,  e.g.,  the 
method  of  inducing  the  separat_on  ("free  interaction").  If  the  final  pres¬ 
sure  is  substantially  higher  than  the  pressure  rise  to  separation,  a  consider¬ 
able  pressure  rise  must  occur  after  separation  over  a  certain  length.  Hence, 
once  the  shock  strength  exceeds  the  magnitude  required  to  provoke  separation, 
the  length  of  the  interaction  region  increases  rapidly. 

CORNER  FLOW 

For  the  case  of  a  compression  comer,  illustrated  in  Figure  k ,  the  effect  of 
separation  is  to  alter  the  flow  geometry  such  that  the  supersonic  flow  will 
be  compressed  in  two  stages,  by  a  compression  shock  at  separation  and  by  a 
reattachment  shock. 
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Figure  4.  Flow  Separation  in  Viscous  Corner  Flow 
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At  reattachment,  the  supersonic  flow  is  compressed  to  its  final  value,  i.e., 
that  value  corresponding  to  the  local  slope  of  the  body  at  that  point.  The 
fluid  with  sufficient  momentum  to  proceed  against  the  reattachment  pressure 
rise  continues  downstream  after  reattachment,  while  that  having  insufficient 
momentum  is  reversed  back  into  the  separated  region. 

EFFECT  OF  STATE  OF  THE  BOUNDARY  LAYER  ON  PRESSURE  PROFILE 


Flow  separations  of  initially  laminar  and  turbulent  boundary  layers  have 
basically  similar  behavior,  but  the  streamwise  scale  of  the  interaction  over 
which  the  pressure  rise  is  spread  is  much  longer  for  laminar  than  for  tur¬ 
bulent  boundary  layers ( Figure  5)-  Therefore  it  is  necessary  to  distinguish 
between  wholly  laminar  boundary  layers  and  wholly  turbulent  boundary  layers. 

TRANSITIONAL  SEPARATION  REGION 

It  has  been  noted  that  the  location  of  transition  relative  to  the  separation 
and  reattachment  points  is  an  important  variable  influencing  the  pressure 
distribution.  Consequently  it  is  of  interest  to  examine  the  effect  of  tran¬ 
sition  in  this  region. 

It  has  been  observed  that  the  transition  point  can  traverse  the  separation 
region  at  almost  constant  Reynolds  number.  (See  Figures  8  and  10,  Reference 
1. )  Analytical  or  empirical  formulations  for  the  determination  of  the  posi¬ 
tion  of  transition  in  the  separated  region  are  therefore  not  available.  One 
can  only  discuss  the  effects  of  transition  on  the  pressure  profile  as  shown 
by  experimental  pressure  distribution. 

The  transitional  separation  region  is  characterized  by  large  pressure  gradi¬ 
ents  near  transition.  This  region  has  been  noted  to  be  unsteady  (Reference 
2).  An  abrupt  pressure  rise  before  the  hinge  line  is  observed  when  transi¬ 
tion  is  near  reattachment  (References  1  to  3)-  This  may  be  observed  in 
Figure  6.  The  plateau  pressure  level  and  the  general  shape  of  the  pressure 
distribution  upstream  of  transition  is  similar  to  those  of  the  laminar  case. 

SOME  THEORETICAL  CONSIDERATIONS 


The  separated  flow  field  can  be  analyzed  by  assuming  a  flow  model  consisting 
of  a  thin,  const ant -pres sure  viscous  mixing  layer  separated  from  a  solid  sur¬ 
face  by  an  enclosed  region  of  low- velocity  air  (dead  air)  (Reference  4).  The 
flow  model  under  consideration  is  represented  essentially  by  the  mixing  of  a 
high-velocity  stream  with  a  low-velocity  stream,  along  the  dividing  stream¬ 
line  between  the  separation  and  reattachment  points  (Figure  7)-  The  dividing 
streamline  divides  the  dead  air  from  the  separated  flow  proper,  and  no  fluid 
enters  or  leaves  the  dead-air  region  (mass  must  be  conserved).  The  uniform 
stream  of  velocity  ue  mixes  with  the  dead  air  and  the  mixing  layer  thickens 
and  grows  parabolically,  beginning  at  the  origin  of  mixing.  The  air  scavenged 
from  the  dead-air  region  is  balanced  by  that  reversed  back  into  the  dead-air 
region  by  the  pressure  rise  in  the  reattachment  zone. 
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Figure  6.  Effect  of  Transition  Location  on  Pressure  Distribution  in  a  Comer  Flow 
(Reference  2) 


SEPARATION  FLOW  PROPER 


In  order  for  a  particle  along  a  streamline  within  the  mixing  layer  to  he  able 
to  overcome  the  pressure  rise  through  the  reattachment  zone  and  pass  down¬ 
stream,  its  total  pressure  P^  must  be  larger  than  the  static  pressure  P1  at 
the  end  of  the  reattachment  zone  (Figure  8,  particle  (a)).  A  particle  such  as 
(b),  however,  has  a  low  velocity  with  a  correspondingly  low  Pt  <  P‘.  Assum¬ 
ing  that  the  boundary-layer  thickness  at  separation  is  zero,  the  asymptotic 
velocity  along  the  dividing  streamline  has  been  calculated  for  two-dimensional 
laminar  flow  to  be  equal  to  u  =  O.587  Ug. 

IMPORTANT  FLOW  SEPARATION  PARAMETERS 

The  main  factors  affecting  the  interaction  between  a  shock  wave  and  a  bound¬ 
ary  layer  are: 

•  Mach  number 

•  Reynolds  number 

•  Shock  strength 

•  Sbace  of  boundary  layer  (laminar,  turbulent  transitional) 

•  Location  of  transition 

For  example,  an  increase  of  Mach  number  increases  the  magnitude  of  the 
plateau  pressure  and  an  increase  of  Reynolds  number  reduces  the  plateau  pres¬ 
sure.  This  can  be  shown  using  the  equation  for  plateau  pressure  for  laminar 
flow 
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which  can  be  rearranged  for  large  Mach  number  to  the  form 
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An  increase  of  the  shock  strength  has  only  slight  effect  on  plateau  pressure 
and  increases  the  upstream  interaction  length  (References  5  and  6). 


The  state  of  the  boundary  layer  influences  the  physical  scale  of  the  inter¬ 
action.  Laminar  separation  extends  over  a  larger  distance  than  turbulent 
separated  flow,  and  the  laminar  pressure  rise  is  more  gradual  as  compared 
with  the  abrupt  pressure  rise  in  turbulent  flow.  It  has  been  recognized 
that  location  of  transition  relative  to  separation  and  reattachment  points 
is  an  important  variable. 
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DISTANCE  ALONG  THE  SURFACE 


Figure  80  Separation  Flow  at  Reattachnient  (Reference  4) 


When  transition  is  close  to  the  separation  point,  the  flow  lacks  two- 
dimensionality  and  is  relatively  unsteady.  Pressure  in  the  transition  region 
frequently  reaches  a  peak  value  which  is  higher  than  in  turbulent  flow,  (in¬ 
cidentally,  the  stability  of  a  separated  boundary  layer  increases  markedly 
with  an  increase  in  Mach  number. ) 

Another  important  parameter  is  the  temperature  of  the  wall.  Wall  cooling  has 
only  a  slight  effect  on  the  magnitude  of  plateau  pressure  but  reduces  the 
interaction  length  (Figure  9)* 

The  effect  of  wall  temperature  on  location  of  separation  point  is  shown  in 
Figure  10  (Reference  7).  The  distance  between  the  beginning  of  interaction 
and  the  separation  point  is  shown  as  a  function  of  the  position  of  the 
beginning  of  interaction,  with  wall  temperature  as  a  parameter.  It  is  noted 
that  as  the  wall  is  cooled,  the  distance  (Xs-Xo)  decreases  and  becomes  zero 
at  some  critical  temperature  ratio.  Simultaneously  the  pressure  at  separa¬ 
tion  becomes  equal  to  the  pressure  at  the  beginning  of  interaction,  and  no 
separation  occurs.  The  limiting  line  on  the  left  indicates  tangency  of  the 
dividing  streamline  to  the  edge.  Also  the  inclination  of  the  dividing  stream¬ 
line  in  the  reattachment  region  is  reduced  by  wall  cooling,  and  as  a  conse¬ 
quence  the  over-all  distance  of  separated  flow  is  decreased. 

THEORETICAL  ANALYSES  OF  SEPARATED  FLOW 


The  complex  and  challenging  problem  of  interaction,  between  a  shock  wave  and 
a  boundary  layer  has  occupied  many  investigators  since  1939,  when  the  inter¬ 
action  phenomena  were  first  observed  by  Ferri  (Reference  8)  during  tests  on 
an  airfoil  in  a  supersonic  tunnel  where  a  favorable  pressure  gradient  was  ex¬ 
pected.  Since  that  time,  numerous  workers  have  investigated  various  facets 
of  the  interference  flow.  Although  our  knowledge  is  still  far  from  complete, 
a  considerable  understanding  of  the  interaction  phenomena  has  been  gained 
over  the  last  two  decades.  However,  in  spite  of  the  persistent  and  long-time 
effort  on  this  subject,  a  satisfactory  and  simple  theoretical  analysis  ade¬ 
quate  for  engineering  applications  does  not  exist.  The  interaction  involves 
the  separation  point,  the  separation  region,  and  re-attachment;  and  each  of 
these  might  be  considered  to  be  one  of  the  unsolved  problems  of  fluid  dynamics. 

An  introductory  discussion  and  a  historical  review  of  the  developments  in  the 
field  of  separated  flow  phenomena  can  be  found  in  References  9,  10  and  11. 

The  simple  case  of  an  oblique  shock  wave  impinging  on  a  flat  plate  has  been 
the  subject  of  numerous  theoretical  and  experimental  investigations.  The 
flat  plate  model  has  the  advantage  of  being  able  to  separate  the  effects  of 
the  interaction  from  those  arising  from  the  pressure  gradient  associated 
with  a  cur/ed  surface.  The  theoretical  study  of  even  such  simple  types  of 
interaction  presents  many  difficulties.  Some  early  investigations  treating 
the  problems  are  described  in  References  12,  13  and  l4. 

One  feature  which  has  received  much  attention  in  theoretical  considerations 
is  the  effect  of  the  shock  wave  on  the  upstream  flow.  It  was  thought  that 
these  effects  could  be  explained  in  terms  of  the  propagation  of  disturbances 
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Figure  10.  Distance  Between  Beginning  of  Interaction  and  Location  of  Separation  as  Influenced 
by  Cooling  (Reference  7) 


upstream  through  the  subsonic  flow  within  the  boundary  layer  near  the  wall. 

A  simple  theory  which  neglected  viscosity  effects  other  than  those  producing 
velocity  profiles  in  the  boundary  layer  was  proposed.  It  was  found,  how¬ 
ever,  that  the  predicted  upstream  effects  were  much  smaller  than  those  observed, 
and  the  differences  between  the  upstream  effects  in  laminar  and  turbulent 
flow  could  not  be  accounted • for  (References  15,  16  and  17 ). 

Another  approach  to  account  for  the  effects  upstream  of  the  interaction  was 
to  consider  the  mechanism  of  the  interaction  between  the  changes  in  thickness 
of  the  boundary  layer  and  the  pressure  changes  which  they  produce.  The 
theories  developed  along  this  line  (References  18  and  19)  seem  to  be  in 
reasonable  agreement  with  experiments  for  weak  shock  waves. 

In  strong  shock  waves,  which  generate  pressure  gradients  in  excess  of  those 
required  to  produce  separation,  the  separation  point  moves  upstream  until 
the  pressure  gradient  drops  to  the  separation  value.  The  problem  of 
boundary- layer  separation  ahead  of  the  impingement  point  is  dealt  with  in 
References  20,  21,  and  22. 

The  mixing  layer  theory  has  been  applied  to  the  interaction  problem  by  Crocco 
(Reference  23).  This  theory  recognizes  a  mixing  between  the  high-velocity 
external  flow  and  the  slower-moving  air  in  the  boundary  layer,  or  the  trans¬ 
port  of  momentum  from  the  outer  stream  to  the  dissipative  stream  as  a  funda¬ 
mental  process  that  determines  the  pressure  rise  that  can  be  supported  by 
the  flow.  Since  the  external  flow  cannot  be  regarded  as  a  known  datum  for  the 
calculation  of  the  dissipative  flow,  it  is  necessary  to  relate  the  thicken¬ 
ing  of  the  boundary  layer  and  the  consequent  deflection  of  the  external  flow 
to  the  pressure  distribution  by  simple-wave  flow  relations.  This  method 
uses  a  generalization  of  the  von  Kannan  momentum  integral  and  treats  the  flow 
as  quasi- one-dimensional  with  properly  defined  mean  velocity  and  temperature; 
it  is  based  on  certain  parameters  characteristic  of  the  boundary- layer  pro¬ 
file.  Interrelationships  among  these  parameters  are  obtained  from  similar 
solutions  of  the  boundary  layer  equations.  The  flux  details  are  lost  by  this 
method,  but  it  is  believed  that  this  simplified  model  preserves  the  main 
features  of  the  interactions  with  the  external  flow. 

Quantitative  disagreement  between  experiments  and  the  theoretical  results 
obtained  by  the  Crocco-Lees  method  for  the  region  upstream  of  separation 
prompted  Click  (Reference  2b)  to  re-examine  this  method.  As  a  result,  an 
improved  correlation  function  based  on  low- speed  theoretical  and  experimental 
data  was  developed.  The  separated  and  re-attached  regions  were  reviewed  and 
a  physical  model  incorporating  the  concept  of  the  dividing  streamline  was  assumed. 
Using  the  Crocco-Lees  approach  and  the  estimated  correction  functions  for 
the  separated  and  re -attachment  regions,  the  problem  of  the  shock-wave  and 
laminar  boundary  layer  interaction  was  analyzed.  The  calculated  pressure 
distribution  showed  a  satisfactory  agreement  with  experiments. 

A  simple  new  method  was  developed  along  the  lines  of  the  well-known  Pohlhau- 
sen  method,  which  is  applicable  to  problems  with  extensive  separation  (Ref¬ 
erence  25).  A  sample  calculation  of  pressure  for  the  laminar  foot  of  the 
interaction  compared  favorably  with  results  obtained  by  the  Crocco-Lees  method. 
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A  more  conventional  integral  method  following  Karman-Pohlhausen  for  the  study 
of  the  interaction  problem  resulting  from  a  swept  planar  shock  wave  and  a 
laminar  boundary  layer  on  a  normal  plane  is  described  in  Reference  26  .  The 
three-dimensional  effects  due  to  sweep  are  resolved  by  considering  the  prob¬ 
lem  in  a  plane  normal  to  the  shock  plane  and  by  considering  the  effect  of  the 
cross- flow  velocity  component  separately.  The  method  is  based  on  a  fourth- 
degree  or  higher  polynomial  for  velocity  and  enthalpy  profiles. 

Another  method  based  on  a  modified  Pohlhausen  approach  with  the  velocity  dis¬ 
tribution  represented  by  a  fifth-degree  polynomial  with  two  undetermined 
parameters  is  presented  in  Reference  27.  The  concept  of  a  dividing  stream¬ 
line  is  used.  This  method  represents  an  improvement  in  pressure  and  shear 
values  over  those  obtained  by  the  one-parameter  method. 

Abbott  et  al . (Reference  28)  describe  a  method  which  is  based  on  solution  of  the 
boundary  layer  equations  by  the  usual  integral  method  with  the  Prandtl- Meyer 
equation  for  the  pressure  rise  at  the  edge  of  the  boundary  layer.  Pre- 
senaration  and  post-separation  regions  are  studied.  Pressure  distribution  is 
cs  culated,  and  the  effect  of  Mach  number,  Reynolds  number,  and  wall  temper? - 
ture  are  established.  However,  the  calculated  temperature  profiles  in  the 
separation  region  are  unrealistic:  the  calculated  flow  temperature  decreases 
as  one  moves  away  from  a  cooled  wall,  thus  the  cool  wall  is  in  the  unlikely 
situation  of  heating  the  separated  region.  This  trend  can  be  ascribed  to  the 
use  of  similar  velocity  and  temperature  profiles  in  both  the  boundary  layer 
and  the  dead-air  region. 

One  of  the  more  recent  methods  is  due  to  Lees  and  Reeves  (Reference  29).  This 
paper  presents  a  theory  which  is  capable  of  including  the  entire  separated 
flow  without  introducing  semi- empirical  feat  ires.  It  is  based  on  the  assump¬ 
tion  that  the  boundary  layer  approximations  are  valid  over  the  entire  viscous 
flow  region.  The  integral  method  with  the  first  moment  of  momentum  in  addi¬ 
tion  to  the  momentum  integral  is  employed.  The  successful  application  of  this 
method  to  separated  and  re-attaching  flows  is  attributed  to  the  proper  choice 
of  the  one- parameter  family  of  velocity  profiles. 

A  numerical  method  of  treating  laminar  separated  flow  was  reported  recently 
by  Nielsen  et  al.  This  method  of  solution  is  based  on  integral  relations  due 
to  Dorodnitsyn.  The  boundary  layer  equations  are  transformed  by.  the  Stewartson 
and  then  by  the  Dorodnitsyn  transformation  such  as  in  Reference  30*  For  axi- 
symmetri.c  flow,  a  Mangier  transformation  is  applied.  The  success  of  this 
method  is  believed  to  depend  primarily  on  the  form  of  the  equation  for  velocity 
profiles  .  The  advantage  of  this  method  over  those  using  a  one-parameter  family 
of  profiles  is  that  it  is  possible  to  obtain  higher  approximations  to  the  solu¬ 
tion  and  thus  permit  assessment  of  how  closely  the  assumed  velocity  profiles 
approximate  the  actual  non-similar  family  of  profiles.  The  solution  of  the 
simultaneous  differential  equations  is  obtained  by  an  IBM  709^  digital  comput¬ 
er.  This  program  which  gives  the  velocity  and  pressure  distribution  through¬ 
out  the  boundary  layer  for  either  two-dimensional  or  axisymmetric  configura¬ 
tions  for  adiabatic  and  non-adiabatic  cases  is  described  in  detail  in 
References  31  end  7* 
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Section  III 


FLOW  PROPERTIES  ALONG  THE  BODY  SURFACE 


LOCAL  PRESSURE  DISTRIBUTION 

The  pressure  distribution  over  the  surface  of  a  vehicle  flying  at  high  Mach 
number  and  low  Reynolds  number  consists  of  inviscid  and  viscous  contributions 
and  is  established  as  a  result  of  the  so-called  viscous-inviscid  interactions. 

Viscous-Inviscid  Interactions 


A  thick  boundary  layer  at  hypersonic  velocity  changes  the  effective  shape  of 
the  body  and  affects  the  inviscid  external  pressure  which,  in  turn,  influences 
the  development  of  the  boundary  layer.  The  resulting  pressure  is  higher  than 
the  inviscid  pressure  in  the  absence  of  viscous  interactions.  A  measure  of 
the  effect  of  the  interaction  on  the  pressure  is  the  interaction  parameter 


X 


\fC 

a 


(XX 


The  viscous  interaction  parameter  X  is  used  for  delineating  the  viscous  inter¬ 
action  phenomena  into  weak  and  strong  interactions  (Ref.  32). 

The  weak  interactions  occur  when  the  deflection  resulting  from  the  growth  of 
the  boundary  layer  is  sufficiently  small.  This  is  usually  expressed  by  the 
relation 

K  =  0^  =  o(l)  or  less 

CU1U  KXO  .  n 

a r  <  eb 

In.  terms  of  the  viscous  interaction  parameter  X,  for  X<1,  the  effect  of  the 
interaction  is  small;  for  X  =  (l)  or  greater,  the  effe-ct  of  the  interaction 
might  be  significant.  Thus,  weak  interactions  appear  at  low  angles  of  attack 
if  the  Reynolds  number  is  high  at  high  Mach  number,  or  if  the  Mach  number  is 
moderately  supersonic  at  low  Reynolds  number.  They  might  also  occur  on  a 
compression  surface  at  sufficiently  high  angles  ol‘  attack. 
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2  dS*  — 

Strong  interactions  occur  when  K  »1  and  >0^  or  if  X»l.  A  strong  pressure 

interaction  region  exists  in  the  flow  over  a  flat  surface  at  low  angles  of 
attack  when  the  local  Mach  number  is  sufficiently  large  and  the  Reynolds  num¬ 
ber  is  sufficiently  small.  The  high  Mach  numbers  and  low  Reynolds  numbers 
which  make  viscous  interaction  important  generally  preclude  turbulence  in 
the  boundary  layer.  The  effects  of  viscosity  on  a  flat  plate  cause  the 
boundary  layer  to  displace  the  inviscid  flow  from  the  plate  by  a  distance 
equal  to  the  boundary  layer  displacement  thickness,  5*.  Since  the  stream- 
wise  pressure  variation  is  dependent  upon  the  local  flow  deflection,  the 
change  of  displacement  thickness  with  distance  along  the  plate  should  be 
found  in  order  to  determine  the  effective  body  shape .  This  new  body  shape 
determines  a  new  pressure  distribution. 

It  is  apparent  that  the  influence  of  any  pressure  gradient  on  the  growth  of 
the  boundary  layer  will  in  turn  influence  the  pressure  distribution.  There¬ 
fore,  a  certain  interdependence  between  the  viscous  and  inviscid  flow  exists 
that  affects  the  establishment  of  an  equilibrium  condition  between  these  two 
flow  regimes . 

Pressure  Distribution  From  Theory 

The  pressure  distribution  of  a  complex  body  flying  at  hypersonic  velocity  is 
obtained  by  numerical  calculations  which  are  rather  lengthy  and  complex.  Some 
approximate  theories  for  determining  the  flow  field  include: 

•  Newtonian  Impact  Theory 

•  Constant  Energy  Solutions 


The  Newtonian  theory  predicts  distribution  for  cases  in  which  the  flow  is 
assumed  to  deflect  parallel  to  the  body.  Constant  energy  solutions  are  based 
upon  the  analogy  of  the  flow  field  generated  by  a  slender  or  blunted  body  and 
a  blast  wave  generated  by  exploding  a  line  charge  along  the  axis  of  symmetry. 
The  blast  wave  solution  is  accurate  far  downstream  and  invalid  locally  near 
the  nose.  It  is  basically  valid  for  a  body  with  zero  thickness. 


A  number  of  expressions  for  the  pressure  distribution  were  obtained  from  theory. 
One  of  the  methods,  as  presented  in  Reference  33?  is  the  basic  tangent  v/edge 
solution  which  leads  to  the  following  expression  for  pressure  distribution  on 
a  compression  surface: 


r 

1 

2 

7(7+1) 

K2  +  7K  Vl  + 

7+1 

K2 
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T 
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where 

K  -  M*e 


(i) 


This  formula  may  be  reduced  for  weak  interactions  and  K<1  to  an  approximate 
form 
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1  +  7K 


P 


and  for  strong  interactions  where  K»1 

P  _  7(7+1)  Tr2 

p  2  K 

1cc 

The  above  formula  is  restricted  to  freestream  Mach  number  greater  than  3, 

j  Cjv 

and  to  a  local  flow  inclination,  0  =  6^  +  of  less  than  20  degrees 

(Ref.  33) 

Bluntness  Effects 

The  solution  for  P/P^,  as  described  above  should  be  increased  linearly  by  a 
bluntness  increment  in  case  the  leading  edge  is  not  sharp.  The  results  of 
blast  wave  analyses  are  utilized  to  calculate  this  contribution.  An  expres¬ 
sion  for  pressure  distribution  due  to  blast  wave  effects  is  given  in  Ref.  34 
as 
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where  Cy  =  0.112  for  air,  Ae  is  the  effective  sweep  angle,  and  Cp  is  the 

drag  coefficient  of  the  leading  edge. 

Equation  (l),  in  conjunction  with  the  bluntness  expression  (2)  was  used  for  the 
predictions  presented  on  Figures  11,  12,  and  13  (for  CD  "  1,2  and  A  =  90°) , 

For  the  particular  experimental  conditions  the  bluntness  contribution  is 
negligible  for  a>10  deg  and  x/d  >1 6  at  Mco<10. 


Boundary  Layer  Displacement  Contribution 

In  the  weak  interaction  regime,  the  boundary  layer  displacement  contribution 
to  K  can  be  found  from  the  classical  compressible  displacement-layer  expres¬ 
sion  presented  in  Reference  35*  For  zero  angle  of  attack: 
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Figure  li.  Pressure  Distribution  Along  Blunted  Flat  Plate,  =  6.96 
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12,  Pressure  Distribution  Along  Blunted  Flat  Plate,  M  -  7,01 


24 


Figure  13.  Pressure  Distribution  Along  Blunted  Flat  Plate,  M  =  14.2 


If  a? iO,  free  stream  quantities  must  be  replaced  by  the  first  estimate  of 
inviscid  local  quantities,  using  for  example  K  =  M^body  and  the  tangent 
wedge  solution,  Equation  1,  or  employing  the  oblique  shock  tables  (Reference 
36).  A  more  complete  treatment  of  this  procedure  is  included  in  subsequent 
discussions. 

Another  way  of  treating  the  case  for  a  ^  0  is  to  regard  the  viscous-induced 
pressure  distribution  as  a  perturbation  on  the  inviscid  distribution  as  follows. 
Once  again,  the  flow  properties  encountered  by  the  boundary  layer  are  approxi¬ 
mated  by  the  inviscid  flow  properties.  The  tangent  wedge  approximation  is 
then  used  to  calculate  the  viscous  contribution,  that  is: 


For  a  more  general  case,  the  viscous  interaction  contribution  as  modified  by 
the  bluntness  included  pressure  gradient  must  be  considered,  fhe  interaction 
regimes  are  classified  in  accordance  with  previous  definitions  of  weak  and 
strong  interaction.  An  expression  for  Kg  is  obtained  following  the  approach 
in  Reference  37.  This  relation  accounts  for  the  boundary  layer  growth  as  in¬ 
fluenced  i»y  the  local  pressure  distribution  proportional  to  a  power  of  dis¬ 
tance,  PjjOCAL  "  xn-  The  following  expression  for  Kg  was  obtained 
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K4  is  determined  from  Figure  14  or  15,  using  the  appropriate  value  for  n. 
n  =  0  (weak  interaction) 

n  =  -1/2  (strong  interaction) 


If  the  inviscid  pressure  gradient  due  to  nose  blunting  is  much  greater  than 
the  viscous  induced  pressure  gradient,  the  applicable  equation  is  (Ref.  37): 
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Figure  15,  Effect  of  Ton  the  Coi 
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•where  corresponds  to  n  =  -2/3  for  blunt  nose. 

Further  investigation  has  resulted  in  the  following  equation  for  the  sharp 
leading  edge  case 
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which  is  valid  for  weak  interaction  with  Xa<l. 

For  strong  interactions  with  XQ,»1,  the  following  equation  presented  in 
Ref.  3T  is  valid 
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where  the  bracketed  expression 


(7) 


If  a  =  0,  freestream  conditions  are  to  be  used.  Expressions  for  the  "G" 

parameter  are: 


G  =  1.7208  (—  +  0.3859)>for  Pr  =  1 

G  =  1.648  +  0.352), for  Pr  =  O.725 

These  equations  are  compared  with  experimental  data  in  Figure  16, taken  from 
Reference  37,  showing  the  appropriate  regions  of  applicability. 

Combined  Effects  of  Incidence  and  Displacement 

An  alternate  procedure  is  to  define  a  total  similarity  parameter  which 
combines  the  angle  of  attack  effect  and  the  flow  displacement  effect  of  the 
ooundary  layer.  This  method  is  valid  if  the  total  deflection  angle  is  less 
than  20  degrees . 
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Figure  16e  Various  Theories  for  Predicting  Boundary-Layer  Displacement  Effects  on 
a  Flat  Plate  in  Air  and  Comparison  with  Experiment  (Reference  37) 
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Figure  160  Various  Theories  for  Predicting  Boundary-Layer  Displacement  Effects  on 

a  Fiat  Plate  in  Air  and  Comparison  with  Experiment  (Reference  37)  (Concluded) 
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where  Kg  is  determined  as  before.  For  Kt>0,  the  tangent  wedge  formula  may 
be  used.  Flat  plates  at  small  negative  angles  of  attack  may  still  show 
positive  pressure  coefficients  on  the  expansion  surface  because  the  vis¬ 
cous  effects  may  be  such  as  to  give  net  positive  flow  deflections.  This  effect 
accounts  for  positive  pressure  coefficients  noted  in  recent  test  data  on  sur¬ 
faces  which  would  normally  be  considered  to  be  in  expansion  flow  (Ref.  38). 

For  such  surfaces,  Newtonian  Theory  would  predict  Cp  =  0  and  a  Prandtl-Meyer 
expansion  would  indicate  Cp<0. 

When  the  expansion  angle  is  greater  than  the  displacement  effect,  Kt<0, 
the  applicable  expression  for  the  hypersonic  similarity  pressure  distribution 
is  (Ref.  39): 
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The  blast  wave  effect  from  leading  edge  bluntness  is  assumed  to  be  additive 
(Figure  ll) . 


From  the  form  of  the  above  equation,  note  that  as  long  as  Kf.  is  negative, 
the  viscous  effect  does  not  dominate.  In  this  case  one  has  a  number  less 
than  one,  raised  to  a  large  power.  The  result  is  a  pressure  ratio 
p 

tJ-  which  is  less  than  unity. 
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Actually  these  results  for  expansion  surfaces  are  less  reliable  than  the  re¬ 
sults  for  compression  surfaces .  No  satisfactory  results  are  available  for 

2 

cases  m  which  K^<  - 

Viscous  effects  were  found  to  be  negligible  on  compression  surfaces  at  large 
angles  of  attack,  especially  at  large  downstream  distances.  That  is,  the 
previous  expressions  for  pressure  distribution  reduce  to: 
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For  large  x/d,  the  blast  wave  effects  diminish  and  P/Pm  approaches  — , 

Pco 

the  oblique -shock  pressure  ratio.  The  oblique  shock  pressure  ratio  may  be 
determined  from  Reference  36. 

Relative  Magnitude  of  Bluntness  Effects  and  Viscous  Effects- 

Most  theoretical  investigations  of  the  interaction  between  the  boundary  layer 
and  the  external  flow  are  based  on  the  assumption  that  the  leading  edge  or 
the  nose  is  infinitely  sharp.  Actual  vehicles,  however,  cannot  tolerate  sharp 
edges  since  some  blunting  must  be  provided  in  order  to  control  the  heat-ing 
rates.  The  relative  magnitude  of  the  inviscid  pressure  generated  by  a  blunt 
forward  Dart  as  compared  with  the  se.i.f-induced  pressure  due  to  viscous  inter¬ 
action  is  discussed  in  the  following  paragraphs. 

Viscous  effects  will  outweigh  any  inviscid  effects  in  the  blunt  region  when 
the  Reynolds  number  based  on  a  dimension  indicating  the  bluntness  of  a  body 
is  of  the  order  of  100.  When  the  Reynolds  number  is  not  too  small,  the 
inviscid  pressure  field  will  dominate  the  forward  part  of  the  body  and  control 
the  growth  of  the  boundary  layer,  and  the  viscous  effects  will  be  negligible . 
Nevertheless,  far  downstream  where  the  inviscid  overpressure  is  small,  viscous 
interaction  might  show  up. 

The  distance  at  which  the  relative  magnitude  of  the  inviscid  pressure  is 
comparable  with  the  pressure  due  to  the  interaction  for  a  body  of  revolution 
with  nose  diameter  d  is  given  by  (Reference  32)* 
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With  an  interaction  parameter  for  strong  interaction  assumed  to  be  and  with 

Cco  =  0(l),  the  self-induced  pressure  becomes  significant  when 
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For  Re^^cM^,  one  can  expect  a  strong  interaction  region  fully  developed  at 
the  nose,  and  an  insignificant  effect  of  the  blunt  nose.  If  Re^  ~  no 
strong  interactions  are  noticeable  at  the  nose.  As  Re^  increases,  the  in¬ 
viscid  effects  spread  downstream  away  from  the  blunt  region,  but  the  inter¬ 
action  might  still  be  important  on  the  slender  part  of  the  body  if  the  inviscid 
overpressure  (P-P^)  is  small,  particularly  because  of  a  rapid  inviscid  pres¬ 
sure  decay  on  axi-symmetric  bodies.  Eventually,  at  high  angles  of  attack, 
the  inviscid  overpressure  is  so  large  that  viscous-induced  pressure  is 
negligible . 
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Table  1  (Reference  4o)  shows  test  flow  conditions  for  blunted  elliptical  cones 
for  a  Mach  number  range  from  5  to  36  •  Calculations  show  that  for  Moo<10,  no 
strong  interaction  can  exist  at  the  nose.  Since  X<1,  only  weak  interaction 
might  be  effective  on  the  slender  part  of  the  body,  and  at  Mm  =  5  these 
effects  are  probably  negligible. 

For  M  >14,  where X»4-,  one  probably  has  strong  interaction  at  the  nose. 

The  expression  to  be  used  in  a  practical  application  for  the  pressure  distri¬ 
bution  calculation  will  depend  on  the  prevalent  flow  conditions.  Another 
parameter  for  estimating  the  relative  magnitude  of  bluntness  viscous  effects 
has  been  established  in  Reference  37.  If 


the  leading  edge  bluntness  effects  are  negligible. 

In  the  case  of  Tcl/lO  the  viscous  contribution  is  small  to  negligible  com¬ 
pared  with  the  blunt,  leading  edge  effects  on  pressure.  For  r<l/2  both  the 
bluntness  effects  and  the  viscous  effects  should  be  considered.  To  compare 
the  different  formulations,  a  number  of  experimental  distributions  from 
Reference  4i  were  used(Figures  11  to  13) .  The  solution  containing  self- 
induced  effects  should  be  used  only  at  »1.  It  is  primarily  useful  on  lee¬ 
ward  surfaces,  where  extensive  viscous  interaction  may  be  expected,  and  for 
high  Mach  numbers  (M^,  >14) .  In  general  the  test  Reynolds  number  was  too  high 
to  induce  substantial  viscous  pressures.  The  inviscid  solution  should  suffice 
under  most  conditions  for  calculating  Po/P^.  The  viscous  pressure  increment 
is  most  noticeable  near  the  nose  where  X  is  larger. and  the  boundary-layer  slope 
is  steeper.  A  calculation  at  the  hingeline  should  be  undertaken  to  determine 
the  significance  of  the  viscous  pressure  contribution.  Figures  11,  12,  and 
13  show  the  decreasing  importance  of  leading-edge  and  viscous  effects  at  large 
downstream  distances  (x/d>l6).  At  these  downstream  positions  the  inviscid 
sharp  wedge  solution  is  the  most  accurate  for  a,>0..  Figure  12  shows  that  the 
inviscid  solution,  P/PTO  =  P  /Pro  +  Pb/Pra,  is  sufficient  for  ct=  +10°  and  -10° 
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LOCAL  FLOW  PROPERTIES 


Actual  vehicle  configurations  can  be  represented  in  many  cases  by  simple 
elements  like  wedges  or  cones . 

Briefly,  it  can  be  said  that  for  two-dimensional  bodies  with  small  bluntness 
and  negligible  viscous  effects,  the  local  pressure  can  be  expressed  as  a  firs4, 
approximation  by  the  pressure  behind  an  oblique  shock  wave  (Reference  36) .  On 
bodies  with  blunted  leading  edge,  pressure  near  the  leading  edge  will  be  higher 
than  that  predicted  by  the  oblique  shock  theory  and  may  asymptotically  approach 
the  oblique  shock  value  far  downstream  from  the  nose. 


For  a  point  relatively  near  a  blunt  leading  edge,  the  shock  loss  theory  of 
Moeckel  (Reference  42)  may  be  used  to  describe  the  properties  of  the  outer 
flow.  The  static  pressure  ratio  Ptt/Pco  must  be  known  beforehand,  whether  from 
tangent  wedge  calculations,  or  blast  wave  theory,  etc.  It  is  assumed  that 
the  streamlines  of  interest  pass  through  the  leading-edge  bow  shock,  and  then 
expand  isentropically  to  the  known  P^/Pq,.  The  leading-edge  shock  plane  has 
its  normal  inclined  to  the  flow  at  an  angle  which  is  approximately  equal  to 
the  sweep  angle  A.  To  find  the  Mach  number  M  corresponding  to  P  ,  one  first 
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notes  the  expression  for  the  ratio  cf  total  to  local  pressure  for  isentropic 
expansion  behind  the  shock  (y=1.4): 
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The  total  pressure  behind  the  shock  and  the  static  pressure  ahead  of  the  shock 
are  related  by: 
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Solving  these  equations  for  M 


M 


o : 


SL 


2 

_  It>  \n  / 

IfcJ  l 


“sl’ 


2  2 
7  m£  cosd  A 


one  obtains: 


c  ..2  2.  /..2  \ 

-  I  6  ripo  cos  A  IMqq  +  5) 

2  2 

M  cos  A  +  5 


-j 


2 


(9) 


T  for  isentropic  flow  is 
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Since  T,  is  constant  for  an  ideal  gas  in  adiabatic  flow,  one  can  calculate  the 
local  ( of )  quantities. 
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Specifically,  the  Reynolds  number  ratio  becomes: 
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In  this  expression  the  perfect  gas  equation,  the  relationship  between  local 
and  total  temperatures,  the  proportionality  of  sound  speed  to  square  root  of 
static  temperature,  and  the  Sutherland  viscosity  law  have  been  employed. 

For  a  point  considered  to  be  relatively  far  from  leading  edge  bluntness, 
oblique  shock  theory  can  be  used  to  obtain  the  properties  of  the  outer  flow 
when  Pa/Poo  is  known  beforehand.  It  is  implied  that  any  contributions  to 
static  pressure  on  the  body,  such  as  from  viscous  interaction  or  even  blunt¬ 
ness  effects,  can  be  represented  as  an  effective  thickening  of  the  body.  Then, 
the  solution  for  Mach  number  behind  the  corresponding  oblique  shock  is: 


As  with  the  shock- loss  method,  all  local  properties  are  now  known.  Reynolds 
number  ratio  may  be  calculated  either  by  equation  10  or  by 
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where  any  appropriate  expression  for  the  viscosity  ratio  may  be  inserted  (for 
instance  the  Sutherland  law) . 


When  bluntness  is  not  zero,  the  actual  will  generally  turn  out  to  be  some¬ 
where  between  M  and  M  .  From  the  pressure  distributions  that  were  exam- 
SL  OS 

ined  in  deriving  the  correlations  given  in  this  report,  it  was  possible  to 
obtain  good  estimates  for  Ma  values  on  surfaces  with  leading  edges  of  various 
bluntness  dimensions.  The  following  formula  for  local  Mach  number  is  suggested 
(Figure  17.) 
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Figure  17.  Biuntness  Effect  on  Local  Mach  Number,  M 
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BOUNDARY  LAYER  THICKNESS 

The  thickness  of  a  boundary  layer  is  an  important  viscous  flow  parameter  fre¬ 
quently  used  as  a  physical  scale  in  investigating  viscous  flow  phenomena. 

The  correlations  of  distance  parameters  describing  the  pressure  distribution 
in  separated  flow  developed  in. the  following  sections  are  expressed  in  terms 
of  the  boundary  layer  thickness.  A  simple  expression  for  boundary  layer 
thickness  was  developed  by  modifying  the  classical  Bias-’ us  laminar  formula: 
the  fluid  properties  used  in  calculating  the  Reynolds  number  are  based  on  the 
reference  temperature,  T*.  For  the  blunt  leading  edge  case,  of  if  large 
viscous-induced  pressure  gradients  are  present,  the  boundary  layer 

/Pw\-% 

thickness  is  reduced  by  the  factor  (=-)  (Reference  43) . 


The  following  equations  for  the  boundary  layer  are  suggested  : 


for  blunt  leading  edge 
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with 


The  predicted  values  of  boundary  layer  thickness  are  compared  in  Figure  18 
with  available  experimental  measurements  from  References  34,  45,  and 
assuming  T  =  T^  for  tests.  It  was  found  that  for  Mc0<6  the  pressure  ratio 

factor  usually  can  be  neglected.  The  boundary  layer  thickness  for  turbulent 
flow  was  expressed  by  a  modification  of  the  classical  formula  from  Reference  47 
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where,  for  turbulent  flow,  in  the  expression  for  77—, 
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Figure  19  presents  a  comparison  of  equation  l6  with  test  data  from  Reference  48. 
Due  to  the  lack  of  experimental  data  no  effects  of  pressure  gradients  for 
turbulent  flow  were  determined. 


REAL  GA3  EFFECTS 

A  blunt  object  flying  through  the  atmosphere  at  hypersonic  speeds  experiences 
severe  heating  from  the  air  which  passes  through  the  standing  bow  shock. 

Part  of  the  kinetic  energy  of  the  flowing  gas  is  transformed  into  thermal 
energy.  This  alters  the  behavior  of  the  air,  which  can  no  longer  be  con¬ 
sidered  a  perfect  medium.  At  high  enough  velocities,  higher  degrees-of- 
freedom  of  gas  species  are  excited.  Vibration,  dissociation,  and  ionization 
of  gas  particles  may  become  important.  The  usual  flow  processes,  coupled 
with  the  rate  processes,  are  responsible  for  these  real  gas  effects. 


ror  increasing  temperature,  significant  quantities  of  energy  are  absorbed  ly 
the  vibrational  mode.  A  molecule  can  aosorb  only  a  certain  maximum  of  energy 
in  this  form  before  its  bond  is  broken  altogether  (dissociation).  Another 
real  gas  phenomenon  is  the  absorption  of  energy  by  electrons.  Complete  de¬ 
tachment  of  an  e.lectron  is  called  ionization. 


From  the  explanation  acove,  it  is  seen  that  "real  gas  effects"  cause  reduc¬ 
tions  in  kinetic  temperatures  by  diverting  energy  into  other  energy  modes. 
This  is  very  important  in  reducing  recovery  temperature.*-  for  high  speed  ve¬ 
hicles  to  tolerable  levels. 
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Figure  18.  Correlation  of  Boundary  Layer  Thickness  for  Laminar  Flow 
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Figure  19. 


Correlation  of  Boundary  Layer  Thickness  for  Turbulent  Flow 


The  adjustment  from  one  equilibrium  state  to  anothei  does  not  occur  instan¬ 
taneously.  An  average  required  number  of  molecular  collisions  can  be  identified 
for  each  energy  mode.  Therefore,  adjustment  times  or  relaxation  times  vary  in 
an  inverse  manner  with  static  pressure.  A  high-speed  ven5.cle  in  the  upper 
atmosphere  may  have  relaxation  times  which  are  considerably  greater  than  the 
time  required  for  a  typical  particle  to  traverse  the  body.  In  this  case,  the 
Mow  is  nearly  "frozen"  rather  than  in  equilibrium,  and  very  little  energy 
exchange  exists,  inis  may  be  treated  as  a  perfect  gas  flow  with  a  constant 
isentropic  exponent  jr responding  to  the  frozen  condition. 

When  the  ratio  of  the  relaxation  time  to  flow  time  is  of  the  order  of  unity, 
non-equilibrium  conditions  might  exist.  If  this  ratio  is  small,  equilibrium 
flow  prevails. 

Considering  compression  by  a  wedge,  frozen  (or  non-equilibrium)  flow  involves 
higher  temperature,  higher  sound  speed,  higher  pressure,  and  lower  Mach  number 
than  equilibrium  flow.  Figures  20  and  21  compare  frozen,  equilibrium,  and 
perfect  gas  pressures  resulting  from  wedge  shocks  for  various  typical  flight 
conditions . 

It  should  bt  noted  that  the  perfecu  gas  and  frozen  gas  assumptions  give  approxi¬ 
mately  the  same  results,  namely  higher  pressures  and  temperatures  than  for 
equilibrium  flow.  The  difference  between  frozen  flow  and  equilibrium  flow 
conditions  increases  with  velocity  and  wedge  angle.  At  an  altitude  of  200,000 
ft  and  a  wedge  angle  8  -  25  deg,  the  real  gas  effects  are  noticeable  at  veloc¬ 
ity  Uq,  — QOOO  ft/sec.  Figure  22  sh cws  the  altitude  effects  on  oblique  shock 
pressure  ratio  for  a  wedge  angle  of  25  degrees.  It  can  be  seen  that  real  gas 
effects  become  important  above  100,000  ft  for  velocities  above  11,000  ft/sec. 

The  equilibrium  values  are  based  on  results  from  Reference  4«. 

For  a  long  enough  vehicle,  one  would  observe  a  streamwise  decay  of  pressure 
from  the  frozen  to  the  equilibrium  value.  Such  a  pressure  distribution  is 
similar  to  the  pattern  that  would  occur  in  the  flow  of  a  perfect  gas  over  a 
slightly  blunted  wedge.  Therefore,  it  is  not  surprising  that  the  shock  shapes 
for  the  two  problems  are  also  similar,  i.e.,  curved. 

For  a  blunt  body,  frozen  or  non-equilibrium  flow  results  in  higher  temperatures. 
The  shock  standoff  distance  s  increased.  Static  pressure,  however,  is  a  little 
below  the  equiliorium  value.  Fig. ire  23  indicates  the  real  gas  flow  regimes  on 
a  velocity-altitude  plot,  for  blunted  surfaces  one  and  ten  feet  in  length. 

It  can  oe  seen  that  flow  conditions  and  corresponding  reaction  rates  for  the 
assumed  gas  model  depend  on  flight  path  and  body  size.  For  a  one-foot  body, 
non-equilibrium  conditions  may  exist  at  a  velocity  of  10,000  ft/sec  above  an 
altitude  of  170,000  ft,  while  for  a  ten-foot  body,  this  altitude  increases  to 
about  210,000  ft.  This  is  because  the  flow  time  on  a  larger  body  is  increased 
and  the  gas  particles  have  more  time  for  adjustment. 

It  would  be  consistent  with  the  work  so  far  on  pre -separation  flow  properties 
to  advise  the  use  of  either  perfect  gas  or  equilibrium  relations  for  each 
problem,  depending  on  whether  perfect  or  equilibrium  flow  is  a  closer 
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Figure  20o  Oblique  Shock  Pressure  Ratio  Pw/Fwvs  5  for  Equilibrium,  Perfect 
Frozen  Gas  Flow 
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Figure  22.  Oblique  Shock  Pressure  Ratio  for  Equilibrium  Flow  and  Perfect  Gas  Flow 
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L  =  10  FT 


approximation.  It  should  be  remembered  that  the  differences  between  perfect 
and  equilibrium  flow  values  decrease  as  deflection  angles  decrease  and  as  the 
altitude  decreases. 

An  accounting  of  real  gas  effects  is  given  in  Reference  50*  Included  therein 
are  some  aspects  of  hypersonic  flow  and  associated  real  gas  effects,  test 
facilities,  some  discussion  of  relaxation  times,  some  results  for  real  gas 
wedge  flows  in  the  atmosphere,  a  description  of  blunt  body  flow,  and  a  rel¬ 
atively  up-to-date  bibliography. 
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Section  IV 


SEMI-EMPIRICAL  CORRELATIONS  FOR  SEPARATED  FLOW 


In  this  section,  the  equations  to  he  used  for  the  prediction  of  separated 
flow  pressure  distributions  are  developed.  The  behavior  of  the  separated 
flow  depends  strongly  on  whether  the  boundary  layer  is  laminar  or  turbulent; 
therefore,  the  discussion  begins  with  an  examination  of  the  transition  phenom¬ 
enon.  A  criterion  for  establishing  the  location  of  transition  is  presented. 
The  next  step  is  to  determine  the  minimum  pressure  rise  required  to  cause  flow 
separation,  and  incipient  separation  criteria  for  both  the  laminar  and  turbu¬ 
lent  cases  are  given.  Finally,  correlation  expressions  are  presented  for  the 
parameters  necessary  for  the  construction  of  the  complete  pressure  distribu¬ 
tions  for  both  laminar  and  turbulent  separation. 

TRANSITION  FROM  LAMINAR  TO  TURBULENT  FLOW 

An  investigation  of  transition  phenomena  was  presented  by  Deem  and  Murphy  in 
Reference  51.  They  collected  large  quantities  of  experimental  data  concern¬ 
ing  transition  of  boundary  layers  in.  high  speed  flows,  and  derived  an  empiri¬ 
cal  expression  for  the  prediction  of  transition  distance.  The  higher  Mach 
number  data  were  taken  mainly  from  the  experiments  of  that  paper ,  and  from  the 
measurements  of  References  52,  53?  and  5^-.  The  scatter  of  the  collected  data 
is  significant;  nevertheless,  there  is  reason  to  believe  that  the  resulting 
expression  is  useful  for  prediction  of  transition.  The  complexity  of  the 
phenomenon  is  shown  by  the  fact  that  the  Deem-Murphy  equation  involves  unit 
Reynolds  number  to  a  fractional  power.  There  is  a  strong  Mach  number  effect, 
with  increasing  Mach  number  tending  to  increase  stability  of  the  boundary 
layer  for  M  >  3.  A  crossflow  factor  must  be  included  when  there  is  sweep. 

Any  effects  of  wall  temperature  were  noted  to  be  surprisingly  slight. 

Certain  comments  concerning  the  final  results  of  Deem  and  Murphy  in  Refer¬ 
ence  51  are  in  order.  That  paper  is  particularly  concerned  with  the  effect 
of  leading  edge  bluntness.  Regimes  of  bluntness  are  defined  and,  in  fact, 
iteration  is  required  to  get  answers .  Furthermore ,  the  expression  for  very 
large  bluntness  includes  a  combination  of  free-stream  and  local  quantities. 

It  is  believed,  however,  that  the  transition  distance  of  the  boundary  layer 
should  depend  only  upon  the  flow  properties  which  it  senses  during  its 
development . 
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Therefore,  the  following  equation  for  transition  distance  is  recommended: 


JL  1 

(feet)  =  ^>38x102+1,94x102  \V:a-3\  2 (cos  A)2 

j(Re/in)a] 


(IT) 


In  this  expression  a  -  subscript  properties  are  the  local  quantities .  These 
local  quantities  are  to  be  calculated  as  well  as  possible  through  the  use  of 
any  available  methods,  such  as  tangent  wedge,  blast  wave,  or  the  Moeckel 
shock  loss  theory  (Reference  42)  when  leading  edge  bluntness  is  comparatively 
large.  (A  more  detailed  discussion  of  local  flow  properties  is  presented 
in  Section  III.)  The  Deem-Murphy  result  and  the  above  result  become  identical 
when  there  is  no  bluntness  and  angle  of  attack  is  zero. 

The  coordinate  XQ  which  identifies  the  beginning  of  the  interaction  is  not 
known  a  priori.  However,  for  vehicles  of  practical  dimensions,  the  percent¬ 
age  difference  between  XQ  and  Xj^r  will  usually  be  rather  small.  Therefore, 
it  can  be  assumed  that  Equation  (17)  predicts  a  turbulent  boundary  layer  if 
Xt  <  X„L  or  Re  <  Re  . 

"D  nil 

To  indicate  the  validity  of  Equation  (l'')  for  prediction  of  the  state  of  the 
boundary  layer  at  the  hinge  line  (or  shock  impingement  point).  Figure  24  has 
been  constructed.  In  this  plot,  the  spectra  of  data  employed  in  devising  the 
interaction  length  and  pressure  level  correlations  for  both  the  laminar  and 
turbulent  cases  have  been  collected  and  compared  with  Equation  (lT)-  The 
Reynolds  number  versus  Mach  number  plot  seems  to  be  appropriate  for  delineat¬ 
ing  the  laminar  and  turbulent  regimes,  but  in  this  case  there  is  no  single 
dividing  line  between  laminar  and  turbulent;  an  X^  (or  alternatively,  an 
Reff/in.)  must  be  specified  because  Equation  (17)  contains  unit  Reynolds  num¬ 
ber.  The  procedure  for  constructing  the  X^  isolines  in  Figure  24  was  to 
select  XHL  and  Re  y  (and  therefore  Re^/in.),  and  solve  for  Ma  from 
aAHL 

Equation  (17)  for  A  =  0  and  Mq,  >  3-  In  this  work  most  interest  is  in  Mq,>3. 
It  is  seen  from  the  plot  that  the  Deem-Murphy  formula  as  suggested  here 
separates  the  believed  laminar  data  from  the  believed  turbulent  data  very 
well . 


The  results  presented  in  Figure  25  indicate  that  for  a  typical  equilibrium 
reentry  glide  trajectory,  executed  by  a  blunted  vehicle,  e .g .  30  feet  in 
length,  control  surfaces  are  likely  to  encounter  both  the  laminar  and  the 
turbulent  boundary  layer.  Very  large  bluntness  was  assumed,  and  the  shock 


loss  theory  was  employed.  Computations  also  show  that  sweep  nss  a  significant 
effect  on  transition  distances.  The  very  high  sweep  of  70°  tended  to  cause 
a  decrease  in  Xt  roughly  25  -  45$. 
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SYMBOL 


XHL  REF  FLOW 

5.5  IN.  83  TURBULENT 
12.0  IN.  84  TURBULEM1 
10.4  IN.  5  TURBULENT 
8.0  IN.  1  TURBULENT 
9.0  IN.  85  TURBULENT 
12.0  IN.  74  LAMINAR 
9.0  IN.  85  LAMINAR 
5.5  IN.  83  TURBULENT 


TEST  DATA 


THE  ATTACHED  FLOW  TRANSITION 
POINT  CONDITIONS  AS  CALCULATED 
BY  THE  DEEM-MURPHY  FORMULA — -* 


10  X  10° 


Re«xt  FORXt  =  XHL  =  5.5IN  = 
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:XHL=  12  IN. 


■  APPROXIMATE  LIMIT  OF  TURBULENT  FLOW  REGIME,  REF.  83 


Figure  24„  Regions  for  Existence  of  the  Laminar  and  Turbulent  Boundary  Layers 


Figure  25«  Transition  Distance  as  a  Function  of  Altitude  and  Angle  of  Attack 


II  should  bo  mentioned  that  the  correlation  formula  presented  here  does  not 
take  account  of  many  factors  which  should  affect  transition  distance  in  some 
manner.  In  particular,,  the  flight  of  a  real  vehicle  through  the  atmosphere 
will  involve  some  conditions  which  are  not  simulated  in  wind  tunnel  tests. 

The  following  factors  should  be  mentioned  relative  to  the  X-^  correlation 
given  above. 

•  Higher  Mach  number  data  is  still  required. 

•  Free  stream  turbuLcnce  Level  is  known  to  have  an  effect  on  ti'ansi- 
t ion  distance  . 

•  Actual  vehicle  configurations  are  more  complex  than  the  wind  tunnel 
flat  plate  models. 

•  It  is  believed  that  there  is  a  stabilizing  effect  of  cooling 
(Reference  l)5) ,  and  a  real  vehicle  may  have  a  protective  cooling 
system . 

•  Surface  roughness  or  protrusions  contribute  to  instability. 

•  R< al  gas  effects,  relaxation  time,  and  flow  time  influence  the  gas 
composition  and  local  flow  properties. 

CONDITIONS  AT  INCIPIFITT  SEPARATION 

It  is  well  established  that  the  existence  of  a  sudden  pressure  rise,  as  at  a 
compression  corner  or  across  an  impinging  shock,  does  not  always  lead  to 
flow  separation.  Gome  minimum  overall  pressure  ratio  for  incipient  separation 
is  required  .  The  first  step  in  the  examination  of  the  interaction  region  is 
the  determination  or  correlation  of  the  quantity: 


There  exists  a  simple  theory  for  the  form  of  the  relation  between  the  pres¬ 
sure  parameters  and  the  other  flow  properties  in  the  separated  region 
(Reference  .  For  this  reason  the  correlations  for  these  pressure  parameters 
are  semi -empirical  rather  than  empirical.  Dynamically,  the  situation  is  that 
of  a  pn-ssnm  rise  overcoming  the  sh°ar  force  between  the  attached  boundary 
layer  and  the  surface.  Since  the  vertical  height  over  which  the  pressure 
differential  acts  is  itself  related  to  the  pressure  increase  through  the 
Prandti-Meyer  relation 
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the  force  balance  is  of  the  form 


CP^U  ~  CP  CpAs]  ~  Cf  A s 


CP  ~ 


For  the  laminar  boundary  layer  C,.  ~  *^Re  Y  ,  and  for  the  turbulent  boundary 

1  aXHL 

layer  a  common  assumption  is  Cf  ~  'W  Re  y  (it  is  considered  that  transi- 

aXHL 

tion  has  occurred  relatively  nearer  X  =  0  than  X  =  X^) .  Therefore  the  data 
for  the  incipient  separation  pressure  are  plotted  as 

C„  f  ReaX,  V*  vs  (MqP-1)  for  laminar  cases,  and 


1_ 

Cp  vs(Mq,2-1)  for  turbulent  cases, 

°INCy  J 

in  Figures  2 6  and  27-  By  establishing  a  straight-line  curve  fit  through  the 
data  in  each  case,  the  following  expressions  result: 


For  laminar  flow: 


CP  (Rea7L  V  =  2*03  2-l)"°*3°6 

Pq,ikc\  ^hl) 


and  for  turbulent  flow: 


■«INC 


=  2.20  (constant) 


Some  experimental  data  shown  in  Figure  26  were  reducted  from  values  for 
separation  point  pressure,  using  a  theoretical  relation  from  Reference  56: 


=  2C 
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Figure  26,  Laminar  incipient  Separation 


Ot/'t 


in 


(m  a  2  -  l) 


Figure  27.  Turbulent  Incipient  Separation 


5^ 


The  theory  in  Reference  56  shows  good  agreement  with  the  experimental  data . 
The  theoretically  derived  expression 


(20) 


where  the  undisturbed- flow  skin  friction  coefficient  is  (Reference  57) 


C  =  .664 
fo 


.45+.5H7  +  -°9  (?~l)M«Pr 
la 


<»>-l 

2 


is  represented  by  a  straight  line  of  nearly  the  same  slope  as  the  correlation 
line,  and  with  only  a  slight  displacement  from  it. 

It  could  not  be  established  that  the  turbulent  flow  follows  a  similar  trend. 
Also,  the  range  of  Reynolds  numbers  in  the  available  data  is  too  small  to 
determine  whether  the  tenth  root  relation  really  is  representative  of  the 
turbulent  data.  Each  expression  does,  however,  correlate  the  available  data 
over  the  ranges  of  test  conditions. 

DEFINITION  OF  THE  INTERACTION  PARAMETERS 


To  describe  the  characteristics  of  the  separated  flow  region  completely  and 
without  ambiguity,  some  definitions  of  the  interaction  parameters  have  to 
be  formulated . 

Figure  28  presents  a  typical  pressure  profile  within  the  interaction  region. 
The  distance  parameters  are  based  on  a  common  reference  line  which  is  defined 
for  various  model  configurations  as  a  geometrical  shock  impingement  point, 
hinge  line  (for  corner  flow),  or  forward  face  of  a  step. 

The  beginning  of  the  pressure  interaction,  x0,  is  defined  as  the  point  where 
pressure  just  begins  to  rise .  The  upstream  spread  of  the  pressure  is 
expressed  by  the  upstream  pressure  interaction  length,  d^,  defined  as  the 
distance  from  the  beginning  of  the  interaction  to  the  reference  line .  The 
free  interaction  length,  is  defined  as  the  distance  from  the  beginning 

of  the  interaction  to  the  point  where  the  pressure  distribution  reaches  the 
straight  section  of  the  pressure  plateau.  The  downstream  interaction  length, 
d2,  is  the  distance  from  the  reference  line  to  the  point  of  intersection  of 
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too  lines  tangent  to  the  pressure  curve,  as  shown,  in  the  downstream  region. 
The  parameter  d^  is  the  distance  from  the  reference  line  to  the  point  of 
intersection  of  two  lines  tangent  to  the  pressure  curve  as  indicated  in 
Figure  28.  The  separation  length,  £QJ  is  the  distance  from  the  separation 
point  to  the  reference  line  (the  separation  point  is  obtained  from  Schlieren 
photographs  as  the  point  of  deflection  of  the  boundary  layer) .  Pressure 
levels  are  defined  corresponding  to  the  previously  mentioned  characteristic 
regions  and  positions.  Therefore,  PQ  and  Pg  correspond  to  pressures  at  x 
and  to  £s  respectively.  Pressure  over  the  separated  region  is  established 
at  a  level  defined  by  Pp.  The  final  downstream  pressure  is  designated  by 

P2- 

SEMI-EMPIRICAL  CORRELATIONS  FOR  LAMINAR  SEPARATED  FLOW 


Flow  separation  occurs  if  the  applied  pressure  rise  is  sufficiently  large, 
i.e.  is  equal  to  or  exceeds  the  pressure  rise  for  incipient  separation. 

Once  this  condition  is  reached  a  characteristic  pressure  distribution 
within  the  separated  region  is  developed. 

In  order  to  approximate  this  pressure  profile,  the  interaction  parameters 
must  be  known .  The  purpose  of  this  section  is  to  present  methods  of  corre¬ 
lating  the  interaction  parameters  into  general  expressions  in  the  form 

Qi  =  f.,  (Mft,  Re, 

These  correlations  are  expressed  in  terms  of  local  properties  as  defined  in 

Section  III,  assuming  that  incipient  separation  conditions  have  been  reached 

or  exceeded.  The  evaluated  experimental  data  cover  a  range  of  Ma  =  1.9  to 

8.4  with  some  data  at  «  14  and  Re  v  =  1.2  x  105  to  1.4  x  10°. 
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ELateau  Pressure  Pp 

One  of  the  most  important  features  of  the  shock  wave  boundary  layer  inter¬ 
action  is  the  plateau  pressure  which  extends  over  the  separated  region.  This 
plateau  pressure  corresponds  to  a  region  of  relatively  constant  slope  of  the 
separated  boundary  layer.  It  has  been  estaolished  from  a  simple  momentum 
equation  argument  (Reference  56)  or  by  using  an  order  of  magnitude  analysis 
(Reference  2),  that  the  plateau  pressure  can  be  expressed  by  the  following 
functional  relation*. 


p^p 


P^  -  P 
P  O 

Y  2 
2  Po 


=  f  (V  Re«X  > 
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or  more  specifically 


(Gp  )p 

La  r 


K(Ma2-l)n 


To  obtain  K  and  n,  the  experimental  values  of  (Cp  )p  ^ReQ,Xoj^-  were  plotted 

versus  (M*2-l)  on  log-log  paper  (Figure  29) .  The  equation  of  the  linear  curve 
fit  through  the  data  is : 


1.56(M02-l)"°*262 


Upstream  Interaction  Length,  d^ 


(21) 


Knowing  the  magnitude  of  the  plateau  pressure,  it  is  still  necessary  to  find 
the  streamwise  spread  of  the  pressure  interaction  region.  Hie  concept  of  the 
free  interaction  states  that  separated  flow  can  be  subdivided  into  a  region 
upstream  of  separation  (which  is  independent  of  the  method  of  inducing  sep¬ 
aration),  and  a  reattachment  region.  One  can  treat  these  regions  separately. 
This  was  the  basis  for  suggesting  separate  distance  parameters  for  upstream 
and  downstream  interaction  lengths .  The  upstream  spread  of  the  pressure  inter¬ 
action  js  expressed  by  the  length  d^. 

It  has  been  observed  that  the  upstream  spread  of  the  -  separated  flow  increases 
with  the  magnitude  of  the  adverse  pressure  gradient  as  expressed  by  the  pres- 
sure  ratio  across  the  shock  wave  or  by  the  magnitude  of  the  flap  angle.  This 
trend  was  indicated  by  experimental  data  for  both  laminar  flow  (Reference  8i.) 
and  turbulent  flow  (Reference  5).  On  the  other  lruid  the  plateau  pressure  In¬ 
creases  slightly  with  the  increasing  flap  angle.  Tims,  U10  magnitude  of  the 
adverse  pressure  gradient  which  has  to  be  negotiated  by  the  boundary  Layer  can 
be  expressed  by  non-dimensional ized  pressure  iifferenco  (P2  -  Pp)/P0,  (P2  -  PQ)/ 
p0,  (pp  “  po)/Po  or  by  parameter  M^.6 .  An  attempt  to  correlate  dj/60  as  a 
function  of  (?2  -  P  )/PpVMa ,  which  was  suggested  in  literature  (Reference  89), 
did  not  lead  to  apsuccessful  and  conclusive  correlation.  A  correlation  of 
di/6G  versus  (P2  -PQ)  /Pq  resulted  in  the  expression 
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ire  29.  Plateau  Pressure  Correlation  for  Laminar  Corner  Flow  (for  Sharp  Leading  Edge  Base  Plate) 


Figure  30  shows  a  plot  of  experimental  data  in  relation  to  the  line  represented 
by  the  above  correlation  equation.  A  correlation  of  d^/60  as  a  function  of 
Mft5  has  led  to  the  following  equation  (Figure  3l) . 


1030 

6o  ‘  ffL-75 


(2.65  )M°6 


The  scatter  of  the  test  points  is  comparable  to  that  of  Figure  30.  Although  the 
experimental  data  inferred  that  the  upstxeam  interaction  length  is  influenced 
by  the  deflection  angle  or  pressure  difference,,  the  evaluated  test  data  did  not 
result  in  an  acceptable  correlation  of  d^/60  as  a  function  c?  these  parameters. 
The  reason  for  this  might  be  the  incompleteness  of  test  data  and  insufficient 
variation  of  evaluated  test  data. 

A  functional  relation  of  the  form 


6 


o 


was  therefore  used  to  correlate  the  experimental  data. 


To  obtain  a  specific  relation,  d]/80  was  plotted  versus  (  7^  -  1  )  (Figure  32) 


\  o  / 

Mach  number  Mawas  noted  for  each  point  and  the  systematic  variation  oi  this 
parameter  across  the  plot  suggested  the  relation: 


where  n  is  the  slope 


of  the  data  for  each  constant  Ma. 
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Figure  30.  Alternate  Upstream  Interaction  Length  Correlation  for  Laminar  Comer 
Flow  (  with  Sham  Leadina  Edae  Base  Plate) 

*  •  W  V  * 


Figure  31 .  Parametric  Piot  for  Alternate  Upstream  Interaction  Length  for  Laminar 
Comer  Flow  (with  Sharp  Leading  Edge  Base  Plate) 


Figure  32.  Parametric  Plot  for  Upstream  Interaction  Length  for  Laminar 
Corner  Flow  (with  Sharp  Leading  Edge  Base  Plate) 
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The  equation  for  F^^)  was  obtained  by  approximating  straight  lines  through 
the  points  of  constant  Mach  number.  When  the  values  of  F(Ma)  were  obtained 
from  the  preceding  equation  and  plotted  versus  Ma  on  log-log  paper  it  was 
found  that  F(Ma)  has  the  form  (Figure  33)  .  From  this  plot  K  and  m 

were  obtained.  Therefore 


(Figure  3*0 


The  resulting  equation  for  d1  which  appears  to  give  best  results ,  is 


9.69  x  lo' 


M 


-4 .  .1 


a 


P 

_P 

P 


(22) 


Free  interaction  Length,  i fi 


The  first  part  of  the  pressure  profile,  from  the  beginning  of  the  interaction 
to  the  plateau  pressure ,  is  called  "free  interaction" ;  its  length  is  defined 
by  the  free  interaction  length,  lf±.  For  the  correlation  of  the  free  inter¬ 
action  length,  the  same  procedure  as  for  the  correlation  of  dp  was  followed 
(Figures  35  and  36) ,  The  resulting  expression  for  the  free  interaction 
length  is 


.47  x  nr  m 


-4.2 


a 


P 

P 

o 


-1 


•49 


(Figure  37)  (23) 
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Figure  36o  Free  Interaction  Length  Correlation  for  Laminar  Corner  Flow  -  Auxiliary  Piot 


e  Interaction  Length  Correlation  for  Laminar  Comer  Flow  (  with  Sharp  Leading  Edge  Base  Plate) 


Downstream  Interaction  Length,  to  Pressure  Rise,  dg 

The  distance  d„  corresponds  to  the  "beginning  of  flow  reattachment  as  defined 
in  Figure  28 .  ^Figure  38  is  -a  detail  sketch  of  this  region . 


Figure  38.  Definition  of  Downstream  Interaction  Parameters 

To  determine  d^,  the  angle  is  first  obtained  from  (f)  =  -vQ,  where  v  is 
the  Prandtl-Meyer  Function.  Table  II  of  Reference  36  yields  v  as  a  function 
of  either  M  or  P/Pt .  The  value  of  vQ  follows  from  and  the  value  of  Vp 
follows  from  Pp/Pt  =  Pp/PQ  x  PQ/Pt .  Once  <f>  is  known,  the  quantity  d^  is 
obtained  from  the  simple  geometrical  relation 

tan  _  L 
tan  dp  +  L 

Therefore : 

d.,  tan  <b  T  d.,  tan  (f> 

T  _  1 _  ,  _  L  _  _ 1 _ _ _  (2k) 

tan  -  tan  <f>  ’  a3  cos  ~  sin  -  cos  8p  tan  <f> 

Downstream  Interaction  Length  to  Peak  Pressure,  &2 

The  location  of  the  final  rise  to  the  pressure  Po  behind  an  oblique  shock  is 
defined  by  the  downstream  interaction  length  dg  (Section  IV).  Two  methods 
have  beer,  used  to  determine  this  parameter. 

Method  1.  The  following  correlation ,  together  with  its  stated  modifications , 
gives  satisfactory  answers  for  any  input  parameters.  Since  dp  is  a  dimension 
of  the  same  order  as  d^,  and  since  Ma&f  indicates  the  strength  of  the  dis¬ 
turbance  cause  by  the  flap  deflection,  plots  of  d^/dp  versus  M^df  were 


TO 


I 


constructed.  The  available  data  had  been  obtained  for  different  flap  lengths, 
and  the  parameter  Cfiap/di  was  introduced  in  order  to  account  for  the  effect 

of  finite  flap  length.  Expressing  the  parameter  'dg/d^  ^fiap/di"  as  a 

function  of  (M^^) ,  the  following  correlation  expression  for  average  values 
is  obtained: 


0.545  -  0.0403  (M^6f) 


(25) 


The  possibility  of  obtaining  unreasonable  answers  for  d2/d^  is  eliminated  by 
imposing  certain  cutoff  points  for  the  above  expression.  These  are  expressed 
mathematically  by: 


For 


For 


dl 


flap  ^ 


0.25: 


d. 


=  .545  - 
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MaSf>  5 
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>  1 


•  25  ,  if  ?  -f--P-  ^  .25 
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.04  (MaSf) 
•  02  (Ma8f) 


<  1 


Method  2.  The  correlation  expression  for  dg  given  below  is  applicable  for  sms 11 

flap  deflection  angles  (6f  <  15  )}  and  flow  conditions  in  the  Mach  number  and 

Reynolds  number  region  M-,  =  1.9  to  8.4  and  Reax  =  1.2  x  10  5  to  1.4  x  10°, 

^  o 


W- 
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respectively.  It  is  again  recognized  that  the  pressure  rise  as  characterised 
by  Sf  should  be  a  significant  factor.  However  8f  itself  is  carried  as  a 
separate  parameter.  To  obtain  a  correlation,  dg/d^  was  plotted  versus 

From  this  a  parametric  dependence  on  8p  was  observed  (Figure  39) . 

A  function  of  the  form 


*'  Uf)  (Ma8f)m 


was  selected  to  fit  the  data .  The  exponent  m  is  the  slope  of  the  parametric 
lines.  The  function  F (8*0  could  then  be  calculated  and  plotted  versus  8f. 
This  resulted  in  F(8^*)  =  K  8pn  (Figure  4o),  where  n  and  K  are  determined  from 
the  graph.  The  final  equation  for  dp  is 


=  8.4  x  10"^ 8  ^  (M  8)1,1''9  (Figure  4l) 

f  w  f 


(26) 


where  8f  is  expressed  in  radians . 

Peak  Flap  Pressure,  P^ 

Because  the  flow  is  turned  more  gradually  when  a  separated  region  exists, 
than  it  would  be  through  a  single  oblique  shock,  the  final  pressure  reached 
on  the  flap  generally  will  be  between  the  isentropic  value  and  the  single 
shock  value  (maximum  entropy  rise).  Therefore,  it  is  suggested  that  for 
laminar  cases  the  average  of  these  two  values  be  used,  with  the  following 
exceptions  (Figure  42):  For  M^>6,  the  single  shock  value  appears  to  be 
better  than  the  average  value  for  the  final  pressui-e  .  For  Ka<6  and  8f  <10°, 
it  is  immaterial  which  of  the  methods  of  pressure  ratio  calculation  is 
employed . 

Separation  Distance,  i 

s 

Some  iterative  numerical  methods  of  analyzing  the  separated  flow  (Reference  7) 
require  the  conditions  at  the  separation  point,  to  start  the  procedure. 


The  actual  separation  point  of  the  boundary  layer  is  defined  by  the  occurrence 
of  zero  velocity  gradient  at  the  wall  in  the  direction  normal  to  the  wall. 


u 


9y 

and  is  correlated  by  the  equation 


=  0  , 


B 


/  1  \  ?.o  . 

-')  KJ 


%  --:.28 


(Figure  43)  (2?) 
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Figure  40e  Downstream  Interaction  Length  Correlation  for  Laminar  Comer  Flow  "  Auxiliary  Plot 


Figure  42„  Peak  Flap  Pressure  for  Isentropic  and  Single  Shock  Compression 
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3.28 


Figure  43.  Separation  length  for  laminar  Flow 
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where 


S 


°B 


•5 .2  x„ 


Separation  Pressure,  P 

s 

The  pressure  level  associated  with  the  separation  position  has  teen  derived 
theoretically  (Reference  56).  It  has  been  found  that  the  pressure  coefficient 
at  separation  can  be  approximated  by 
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where  is  the  skin  friction  coefficient  for  the  undisturbed  region.  Using 
the  expression  for  CpQ  from  Reference  30>  the  following  relation  was  obtained. 


r  t  1 1 

co-1 

f 1.328 

0.45  +  0.55  Tjr^-  +  0.09(7-l)M2a  Pr  2 

2 

R 


1 

2 


(««  -1) 


ax 


o 


(28) 


The  pressure  at  the  separation  point  for  laminar  flow  has  been  correlated  for 
various  configurations  in  Reference  6,  resulting  in  a  second  expression  for 


0.91  (m2  J-0-28t 

K)* 


(29) 


A  comparison  between  the  two  expressions  may  be  seen  in  Figure  44.  The  good 
agreement  between  these  two  expressions  tends  to  indicate  the  validity  of  the 
theoretical  relation . 


The  Effect  of  Finite  Flap  Span 

Pressure  distribution  in  the  separated  region  over  the  flap  is  influenced  by 
finite  flap  span.  Some  data  for  partial  span  flap  pressure  distributions  are 
presented  in  References  58,  59  and  60.  Figure  45  shows  the  configuration 
considered . 
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It  is  desired  to  construct  some  empirical  correlation  for  an  average  flap 
pressure  defined  by  the  equation: 


y 

/  (pg  _p)  dy  =  cpg  yf  =  (pg  "  pav^f 


p  =  P.  (1  -  c) 

av  £ 


The  pertinent  span  parameter  -was  not  varied  in  the  available  data.  The 
average  value  of  C,  which  indicates  average  percentage  loss  in  flap  force 
for  the  evaluated  tests,  was  1 4$>.  Attempts  to  identify  Mach  number,  Reynolds 
number,  or  position  effects  on  C  values  were  unsuccessful ,  The  flap  aspect 
ratio  of  the  test  model  was 


2 

Yf 

AR  =  . -  =  1.33 

yf  flap 

Presumably,  C — >0  as  AR — >  oo,  and  percentage  flap  force  loss  C  should 
increase  monotonically  to  unity  as  AR  decreases,  i.e.,  as  the  deviation  from 
two-dimensionality  increases .  In  the  absence  of  more  information  for  partial 
span  flap  effect,  an  approximation  for  (l-C)  for  design  purposes  is  indicated 
by  the  broken  line  on  Figure  46. 

SEMI-EMPIRICAL  CORRELATIONS  FOR  TURBULENT  SEPARATED  FLOW 


The  distance  and  pressure  parameters  as  shown  on  Figure  28  are  again  used  to 
describe  the  pressure  profile  .  The  methods  used  to  obtain  the  correlations 
for  turbulent  flow  are  virtually  identical  to  those  used  in  the  laminar 
problem.  The  correlations  cover  a  Mach  number  and  Reynolds  number  range  of 
Ma  =  1.5  to  6.2  and  Re  =  1  x  10 ^  to  6  x  10&,  respectively. 


Plateau  Pressure,  PT 


The  plateau  pressure  is  defined  as  the  approximately  constant  pressure  reached 
in  the  separated  flow  region .  Using  the  parameters  suggested  by  theory  and 
by  previous  correlations,  (Cp  )p  was  taken  to  be  a  function  of  Reax  and  Ma. 

^  w  o  2 

The  correlation  was  obtained  by  plotting  (Cp^p  Re^  versus  (M^  -  l) . 

This  led  to  an  expression  for  (^p^p  the  form 


) 

P 


Kofi*  -l)n 

KT 
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where  K  and  n  were  determined  from  the  plot  in  Figure  4?  "by  finding  an 
equation  for  a  straight  line  drawn  through  the  data.  The  following  expres¬ 
sion  for  (Cp  )p  was  obtained. 


^CPa  ^ 


1.91  (*4  -i)"‘3°9 


y  m' 
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(Re  ) 
'  ax 
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Upstream  Distance  Parameter,  d^ 


The  quantity  d-^  was  correlated  by  assuming  the  functional  relation: 


/K,  , 


P  -  P 
P  o 


(30) 


The  distance  d1  was  non-dimensionalized  by  the  boundary  layer  thickness,  80, 
where : 

.15^x0 

8^  = 


'  * 

Re  . 
ax 

o  . 
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From  a  plot  of  (d]_/fi0)  versus  (Pp/PQ-l)  a  function  of  the  form 

cL  /  P  \  n 

_  .  F(Ma)(-^  -1) 


was  inferred,  and  both  n  and  F'(Ma)  were  obtained.  The  reduced  test  results 
are  shown  in  Figure  48.  The  resulting  final  equation  for  d]_  has  the  form 


=  1.1  x  106 


M 


■  'U  - ) 


8.55 


(31) 
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Figure  48.  Upstream  Interaction  Length  Correlation  for  Turbulent  Corner  Flow 
(with  Sharp  Leading  Edge  Base  Plate) 


Free  Interaction  Length, 

The  determination  of  ^fi»  the  free  interaction  length,  followed  a  procedure 
identical  with  the  determination  of  dj_.  A  functional  relation  of  the  form 


n 


where  F(Ma)  =  KMqP,  was  found .  The  exponent  and  constants  were  obtained  from 
Figure  49* 

The  resulting  equation  for  is 


1.84  x  10 


4 


8.40 


(32) 


Downstream  Interaction  Length  to  Peak  Pressure,  d2 


The  downstream  distance  parameter  dg  indicates  the  location  where  pressure 
reaches  the  final  peak  flap  value,  as  indicated  in  Figure  28. 

Method  1.  A  functional  relationship  in  the  form 


■  F^f)  m 

was  used  for  correlation.  A  plot  of  d2/d^  versus  Ma  led  to  a  dependence  on 
flap  angle 

F(8f)  =  K  (8ff 

The  final  relation  was  therefore  of  the  form 


K(8f)n  (MaSf)m 
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Figure  49.  Free  Interaction  Length  Correlation  for  Turbulent  Comer  Flow  (with 
Sharp  Leading  Edge  Base  Plate) 


The  best  fit  to  the  data  was  given  by 


~~  =  2. a  x  io"5 

di 


M 


a 


1  4.934 


(Figure 


where  8f 
8f  <  15% 


is  expressed  in  radians .  The  formula  is  based  on  test  results  for 
=  1.5  to  6.2  and  =  T  x  10^  to  5-6  x  10^. 


Method  2.  An  alternate  approach  leads  to  an  expression  for  d  /cL  as  a 
function  of  and  flap  length  ratio,  of  the  form  1 


d2/di 


'flap 


=  1.16  -  0.33  Ma8f 


(34) 


For  meaningful  results  the  following  limits  are  recommended: 
For  C, 


flap 


>  l 


d2/d±  =  1.16  -  0.33  Ma  8f. 


For 


'flap 

d. 


<  0.25:  d2/d  =  0.58  -  0.165  Ma  8 


For 


where 


Ma  8f  >  2.4:  d2/d1 


r 


1,  for  £  1 

dl 


=  0.37, 


r  = 


flap 


,  for  1  £ 


'flap 


u  .d ; 


0.25,  for 


'flap  < 


s  0 .25 


This  formulation  may  be  used  for  a  wider  range  of  flow  conditions  than  that 
of  Method  1. 
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Figure  50,.  Downstream  Interaction  Length  for  Tui 
(with  Sharp  Leading  Edge  Base  Plate) 
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Comer  Flow 


Downstream  Interaction  Length  to  Pressure  Rise,  dg 


Examination  of  the  test  data  used  in  this  study  led  to  the  conclusion  that  dg 
can  he  taken  as  zero  for  the  turbulent  case .  It  will  be  recalled  that  this 
parameter  indicates  the  beginning  of  the  rise  from  the  plateau  level  to  the 
final  flap  pressure . 


Peak  Flap  Pressure,  Pg 

The  less  extensive  separation  region  typical  of  turbulent  flow  allows  the  use 
of  the  single  shock  method  to  determine  the  pressure  level  Pg.  This  pressure 
value  can  therefore  be  obtained  from  Reference  36  using  the  local  Mach  number 
value  Ma  as  suggested  in  Section  III. 


Conditions  at  Separation,  Ps  and  is 

0 

The  separation  position,  Xs,  again  defined  by  g— 


w 


=  0,  may  be  approxi¬ 


mated  by  d]i  for  turbulent  flow.  The  associated  pressure  coefficient  was 
found  to  be  (Figure  k$) 


<Cp  >s 
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(Reference  6) 


(35) 
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Section  V 


FLOW  SEPARATION  ON  VARIOUS  CONFIGURATIONS 


The  semi- empirical  correlations  for  separated  flow  described  in  the  previous 
section  were  derived  from  experimental  results  for  two-dimensional  bodies. 
Although  two-dimensional  bodies  as  vehicle  elements  are  of  considerable  im¬ 
portance,  actual  vehicles  are  more  complex  and  the  flow  is  three-dimensional 
in  character.  In  this  section  separated  flows  on  fin-plates,  flat  surfaces 
with  compression  flaps  and  endpiates,  and  delta  wings  and  their  variations 
are  discussed.  Three-dimensional  features  of  flow  about  pseudo-two¬ 
dimensional  bodies  of  finite  size  are  discussed  qualitatively  with  the  aid 
of  visual  data. 

FIN-PLATE  CONFIGURATION 


Vertical  fins  and  control  surfaces  are  employed  on  various  types  of  hyper¬ 
sonic  vehicles,  such  as  those  designed  for  aerodynamically  maneuverable 
reentry  missions.  Shock  waves  and  pressure  rises  induced  by  such  devices  at 
finite  incidence  will  be  transmitted  to  the  adjacent  surfaces  on  the  vehicle 
thus  modifying  the  flow  fields  there.  The  problem  of  a  fin-generated  shock 
sheet  impinging  normally  upon  a  surface  with  a  boundary  layer  is  of  interest 
from  both  an  aerodynamic  and  thermodynamic  point  of  view.  The  abrupt  pressure 
rise  across  the  shock  in  real  flow  is  diffused  and  spread  over  some  finite 
distance.  Figure  52  shows  a  simplified  model  of  the. interaction  region. 

Choice  of  Coordinate  System 


In  attempting  to  devise  a  set  of  correlations  for  the  fin-plate  pressure  dis¬ 
tribution,  as  was  done  in  the  deflected  flap  case,  it  must  be  realized  that 
one  is  dealing  with  a  basically  three-dimensional  flow.  The  impinging  shock 
sheet  is  considered  to  be  perpendicular  to  the  plate  surface,  but  it  is  not, 
in  general,  normal  to  the  local  velocity  vector.  The  question  arises  as  to 
whether  for  the  correlation  purpose  one  should  consider  pressure  distribution 
in  the  initially  streamwise  direction  or  proceed  in  a  direction  normal  to  the 
impinging  shock. 

Even  though  the  flow  behind  the  shock  is  not  aligned  in  the  original  stream 
direction,  the  one -dimensional  assumption  may  be  a  good  approximation  (as  long 
as  the  flow  deflection  angle  is  not  too  large) . 
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IMPINGING  SHOCK 


Figure  52.  Detail  of  Shock  Impingement  Interaction 


Correlations  in  the  direction  normal  to  the  shock,  assuming  that  the 
tangential  velocity  component  does  not  change  across  the  shock  and  does  not 
have  significance  for  the  external  inviscid  flow  (Reference  26) ,  may  not  he 
valid  for  the  boundary  layer,  where  the  tangential  velocity  varies. 

Analysis  of  Experimental  Results 

Experimental  Procedure.  The  model  employed  in  the  experimental  study  of 
fin -pi  ate  interaction  (Reference  87)  is  shown  in  Figure  53»  The  fin  or  wedge 
which  generates  the  impinging  shock  has  a  semivertex  angle  of  15°.  Actually, 
data  have  been  taken  on  models  of  four  different  configurations ,  as  there  are 
tall  and  short,  sharp  and  blunted  fins.  Pressure  taps  are  scattered  on  one 
side  of  the  plate  surface  and  one  side  of  the  fin  (the  other  half  of  the  model 
is  instrumented  with  thermocouples) .  Unit  Reynolds  numbers  are  varied  moder¬ 
ately,  with  the  Deem-Murphy  criterion  indicating  that  boundary  layers  are 
probably  laminar  (Section  4).  Considerable  variation  of  the*  local' Mach 
number  is  accomplished  by  angle  of  attack  variation.  Because  of  the  pressure 
instrumentation  arrangement,  the  most  informative  runs  were  those  of  Mro  =  5 
and  a  =  5%  i.e.,  Mq,  =  4.5. 

Shock  Wave  and  Boundary  Layer  Effects.  Figure  53  depicts  part  of  the  system 
of  shocks  that  is  produced  when  the  model  1.  pitched  at  some  compression 
angle  of  attack  in  the  high  speed  flow.  One  shock  sheet  extends  downstream 
from  the  leading  edge  of  the  plate;  it  will  intersect  another  shock  sheet 
generated  by  the  fin,  assuming  the  fin  is  tall  enough.  The  intersection 
ideally  results  in  the  production  of  a  new  shock  system  which  in  general 
undergoes  subsequent  impingements  and  intersections . 


It  appears  that  the  simple  flow  regions  as  noted  on  Figure  53  are  adequate  as 
initial  and  final  conditions  for  the  interaction  problem  examined  here.  In 
particular,  the  "a"  region  is  considered  to  be  the  result  of  turning  the  free 
stream  to  the  plate  surface  inclination.  The  "w"  region  is  considered  to  be 
the  result  of  turning  the  free  stream  to  the  direction  of  the  fin  side.  The 
"2"  region  is  considered  to  be  the  result  of  turning  the  "oc"  flow  to  the  di¬ 
rection  of  the  fin  side.  Therefore  the  initial  and  final  conditions  for  the 
study  of  the  plate  pressure  distribution  through  the  impinging  shock  are  the 
conditions  of  region  "a"  and  the  conditions  of  region  "2"  (Figure  53).-  This 
seems  to  work  fairly  well  in  most  cases.  There  is  generally  some  pressure 
overshoot  before  the  region  "2"  conditions  are  reached.  Also,  the  viscous  in¬ 
duced  deviation  of  the  pressure  ratio  Fa/?co  *r0D1  unity  in  zero  angle  of  at¬ 
tack  cases  is  best  represented  by  isentropic  compression  for  the  Mach  number 
range  tested. 

It  will  be  seen  that  the  boundary  layer  acts  at  all  times  to  smooth  out  pres¬ 
sure  differences.  Also,  the  boundary  layer  on  the  fin  complicates  the  ex¬ 
amination  of  the  interaction  on  the  plate;  an  effective  curvature  of  the  fin 
surface  introduced  by  the  parabolic  (8  ~  \[x.)  growth  of  the  boundary  layer 
causes  the  shock  to  curve. 
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a  ~  FLOW  COMPRESSED  BY  LEADING  EDGE  SHOCK 
w  -  FLOW  COMPRESSED  BY  FIN  SHOCK 

2  -  FLOW  COMPRESoED  BY  LEADING  EDGE  SHOCK  AND  THEN 
BY  FIN  SHOCK 


Figure  53.  Fin-Plate  Configuration  with  Intersecting  Shocb  and  Simple  Flow  Regions 
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Plate  Pressure  Distributions  for  Sharp  Fin.  A  limiting  factor  in  the  analysis 
of  the  tests  is  the  spacing  between  pressure  taps;  more  taps  are  required  in 
order  to  get  a  detailed  picture  of  the  pressure  distribution  through  the 
region  of  the  impinging  fin -induced  shock.. 


A  combined  or  aggregate  pressure  distribution,  for  deriving  the  maximum  of 
information  from  the  data,  can  be  constructed  as  follows.  The  position  of 
the  impinging,  fin -produced  shock  is  estimated  as  well  as  possible,  as  the 
locus  of  the  straight  oblique  shock  that  would  be  produced  by  turning  the 
"a"  flow  through  the  angle  5  .  Then,  all  the  pressure  tap  readings  are 
plotted  against  streamwise  or  normal -to -shock  tap  distance  from  the  assumed 
shock  locus,  on  a  single  set  of  axes. 

Figures  5^  and  55  indicate  the  effect  of  fin  side  boundary  layer.  The  same 
pressures  are  plotted  in  Figure  55  as  in  Figure  5^,-hut  this  time  against 
slighfly  different  relative  streamwise  coordinates,  which  were  obtained  by 
constructing  a  shock  locus  modified  by  introducing  slight  shock  curvature 
expressed  by  Xs(y) ,  which  gives  each  pressure  tap  a  new  relative  position. 
Figure  56  uses  the  normal -to-shock  direction  for  constructing  the  aggregate 
plate  pressure  distribution.  It  should  be  compared  with  Figure  57,  which 
shows  the  distribution  in  the  streamwise  direction:  the  pressure  profile 
normal  to  the  shock  shows  less  scatter  and  its  scale  is  shorter.  Although 
the  plate  pressure  plots  indicate  that  there  is  a  region  characterized  by 
intermediate  pressure  values ,  it  is  not  possible  to  identify  a  plateau  level . 

Sharp  Fin  Pressure  Distribution.  Figures  5ft  and  57,  which  show  plate  pressure 
distributions  for  the  small  and  large  fins,  respectively,  indicate  that  plate 
pressure  is  not  influenced  by  fin  height.  Figure  53  indicates  that  pressure 
distributions  up  the  sides  of  the  fins  will  show  variation  between  the  ,,2,, 
values  and  the  "w"  values.  Representative  results,  shown  in  Figures  58  and 
59,  indicate  that  the  pressure  rise  from  ?w  to  P 2  through  the  shock  inter¬ 
action  is  spread  out  over  an  appreciable  length.  Figure  60  shows  that  the 
pressure  distribution  on  the  fin  wall,  for  a  =  0,  has  a  definite  drop  toward 
the  corner.  The  Pp  value  based  on  displacement  effects  is  not  reached  at  all. 

No  conclusions  are  here  stated  concerning  the  corner  flow  between  the  fin  and 
the  plate,  as  very  dense  instrumentation  is  required  for  a  study  of  this 
region.  (See  Reference  6l.) 


HLuntness  Effects.  The  effects  of  fin  leading  edge  bluntness  are  examined 
through  the  plots  of  Figures  6l  and  62.  Figure  59  shows  that  the  use  of  the, 
ideal,  sharp-leading-edge  oblique-shock  locus  is  no  longer  adequate,  for  it 
results  in  a  greatly  distorted  representation  of  the  plate  pressure  dis¬ 
tribution.  Figure  52  indicates  higher  pressure  close  to  the  blunted  leading 
edge  as  a  result  of  increased  strength  of  the  detached  bow  shock.  The  invis- 
cid  pressure  Pg  is  exceeded  at  the  station  closest  to  the  nose. 


Remarks  on  a  High  Mach  Number  Run.  In  Figures  63  and  60,  pressure  distribu¬ 
tions  for  the  model  at  zero  angle  of  attack  and  =  8  are  presented. 

It  is  seen  that ,  because  of  the  viscous  interaction  effect ,  the  pressure  on 
the  plate  upstream  of  the  fin  is  substantially  larger  than  the  free-stream 
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Figure  54»  Aggregate  Pressure  Measurements  in  Free  Stream  Direction  Along  Plate  with  Impinging  Norniai  Shock 
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Figure  57.  Aggregate  Pressure  Measurements  in  Free  Stream  Direction  Along  Plate  with  Impinging  Normal  Shock 
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Figure  53.  Pressure  Distributions  Along  the  Fin  Side  (Small  Fin) 
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Figure;  59.  Pressure  Distributions  Along  the  Fin  Side  (Large  Fin) 
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Figure  60.  Pressure  Distributions  Along  the  Fin  Side  (High  Mach  Number) 
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Figure  61 .  Aggregate  Pressure  Measurements  in  Free  Stream  Direction  Along  Plate  with  Impinging  Normal 
Shock  (Blunt  Fin) 


Figure  62.  Pressure  Distributions  Along  the  Fin  Side  {Blunt  Fin) 
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value.  However,  the  dominant  pressure  downstream  is  the  pressure  that  would 
be  reached  by  a  single  shock  turning  from  the  free  stream  to  the  direction 
of  the  fin  wall,  i.e.,  Pw  rather  than  P2.  Apparently  the  viscous  displace¬ 
ment  compression  is  insignificant  for  these  configurations. 

The  high  local  Mach  numbers  in  this  case  mean  that  shock  surfaces  are  in¬ 
clined  rather  acutely  with  the  original  free  stream  direction,  and  a  longer- 
model  with  more  downstream  measuring  points  would  be  required  to  obtain  a 
good  view  of  the  shock  impingement  interaction. 

Cor.clusions.  Figures  54  to  60  indicate  that  the  interaction  between  the 
plate  and  fin  gives  neither  a  uniform  pressure  rise  nor  a  single  intermediate 
plateau  pressure.  Projection  of  measurements  describing  the  interaction  re¬ 
gion  on  a  plane  normal  to  the  shock  plane  seems  to  be  feasible.  However,  a 
more  complete  set  of  measurements  on  a  heavily  instrumented  model  would  be  re¬ 
quired.  Understanding  of  the  impinging  shock  interaction  is  necessary  to  al¬ 
low  correlation  of  the  average  pressure  and  overall  spread  of  the  interaction. 

FLOW  ABOUT  A  FLAT  PLATE 


Simple  shapes  are  considered  two-dimensional  if  their  span-to-chord  ratio  is 
very  large.  However,  actual  vehicles  are  of  finite  dimensions,  and  body  ends 
cause  deviation  of  an  actual  flow  from  the  idealized  two-dimensional  flow. 

Flat  Plate  without  End  Plates 


The  tips  of  a  finite  "two-dimensional"  body  in  a  subsonic  stream  influence 
the  pressure  distribution  on  the  entire  body.  At  supersonic  speed  finite 
span  effects  are  also  important,  but  in  a  different  way.  For  a  flat  plate  at 
positive  angle  of  attack,  higher  pressure  on  the  compression  side  tends  to 
equalize  with  the  lower  pressure  on  the  top  surface.  However,  this  pressure 
equalization  and  propagation  of  disturbances  is  confined  to  the  Mach  cone. 

The  resulting  flow  around  the  longitudinal  edges  with  pressure  differential 
leads  to  the  generation  of  vortices,  similar  to  those  existing  in  subsonic 
flow,  which  are  shed  downstream  and  stay  within  the  Mach  cone.  These  vortices 
emanating  at  the  tips  of  the  body  generate  a  typical  circulatory  flow  indi¬ 
cated  by  steamlines  in  Figure  64. 

A  simple  and  illustrative  example  of  end  effects  on  a  two-dimensional  flow  is 
a  flow  about  a  flat  plate.  Figures  65  to  71  are  taken  from  the  data  of  Ref¬ 
erence  74.  Figure  65  shows  an  oil  flow  pattern  for  a  flat  plate  (with  a 
sharp  wedge  as  the  leading  edge)  at  zero  angle  of  attack.  This  figure  indi¬ 
cates  a  distinct  flow  dividing  line  radiating  inboard  from  the  tip,  with  the 
flow  outside  of  this  line  directed  toward  the  side  edge,  and  the  flow  inboard 
of  this  line  directed  toward  the  centerline.  Near  the  centerline  the  flow 
direction  approaches  the  free  stream  direction.  As  a  result  of  higher  pressure 
on  the  bottom  surface,  a  coiled  vortex  is  generated  at  the  tip  and  the  flow 
spearates  near  the  side  edge.  The  vortex  flow  re-attaches  forming  an  attach¬ 
ment  line.  Along  this  line  the  flow  divides,  turning  either  inward  or 
outward,  (Figure  66). 
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Figure  65.  Flat  Plate  Oil  Flow  Study  (End  Plates  Off), 


Figure  67  shows  a  fldw  pattern  on  the  compression  surface  of  a  flat  plate  at 
an  angle  of  attack  of  25°  .  High  pressure  on  the  compression  surface  results 
in  a  conical  flow  directed  outboard.  Figure  68  shows  .a  flow  pattern  for  a 
ramp  at  20°  positive  angle.  Flow  on  the  forward  plate  at  zero  angle  of  attack 
is  basically  the  same  as  that  shown  in  Figure  67.  Downstream  of  the  hinge 
line,  slight  deflection  of  flow  pattern  by  compression  on  the  ramp  is  visible. 
Figure  69  shows  a  plate  inclined  at  25°  with  a  ramp  at  30°  Conical  flow  on 

the  forward  plate  is  noticeable.  An  irregular  separation  line  and  strong 
three-dimensional  effects  due  to  high  flap  angle  are  seen  upstream  of  the 
hinge  line. 

Flat  Plate  with  End  ELates 


The  effect  of  side  plates  on  flat  plate  flow  is  shown  in  Figure  70.  The  re¬ 
gion  near  the  side  edges  indicates  the  boundary  layer  displacement  effect  and 
possibly  the  effect  of  an  oblique  shock  generated  by  the  siu.t=  plate.  The 
surface  flow  lines  indicate  that  this  is  a  region  of  high  shear  stresses  and 
high  pressure. 

Figure  71  shows  a  flat  plate  with  side  plates  at  zero  angle  of  attack  and  a 
flap  at  an  angle  of  10°.  The  forward  plate  flow  is  essentially  similar  to 
the  flow  of  Figure  70 ;  however,  the  flap  region  shows  a  rather  complex  flow 
with  three-dimensional  features  and  possibly  a  multi-vortex  system. 

The  separated  flow  beneath  the  boundary  layer  near  tb<  wall  is  subsonic,  and 
therefore  allows  the  propagation  of  disturbances  in  all  directions.  Effects 
due  to  end  conditions  and  the  effects  of  finite  plate  span  are  expected  to 
be  noted  in  both  the  viscous  and  inviscid  flows.  The  resulting  flow  field  is 
very  complex.  An  experimental  comparison  between  pressure  distribution  over 
a  flapped  flat  plate  with  and  without  end  plates  for  identical  conditions  has 
been  undertaken  (References  59  and.  60) •  Pressure  measurements  we  e  taken  in 
both  the  longitudinal  and  transverse  directions . 


For  =  5 ,  a  -  0,  =  10°  ,  20°  and  Re^/ft  =  1.1  x  10  ,  and  6.6  x  10 

(References  59  and  60) ,  a  number  of  general  trends  may  be  noted  for  the  end 
plates  "on"  and  "off"  configurations.  The  streamwise  pressure  level  on  the 
flap  at  a  flap  deflection  angle  6f  of  10°  is  lower  with  the  end  plates  on  for 
Rejft  -  1  .1  x  10°  and  3.3  x  10°,  At  Re^ft  =  6.6  x  10°  this  trend  re¬ 
verses  and  die  pressure  level  on  the  flap  for  the  end  plates  "off"  condition 
is  lower.  The  overall  differences  are  seen  to  be  small  (Figure  72). 

For  the  same  Mach  number  and  range  .  Reynolds  numbers ,  identical  trends  are 
observed  for  the  6 %  -  20°  cases  (Figure  73).  If  the  plateau  pressure  region 

extends  over  a  significant  distance,  the  effects  of  end  plates  may  become 
important.  For  the  previous  flow  conditions  at  6^  =  20°  and  10°  the  end 
plates  "on"  plateau  pressure  is  slightly  higher  for  all  but  the  Re^/ft  = 

6.6  x  10°  cases.  For  the  higher  Reynolds  number  the  difference  in  the  pres¬ 
sure  level  between  the  end  plates  "on"  and  "off"  condition  disappears 
(Figures  72  and  73). 
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Figure  67.  Flat  Plate  Oil  Flow  Study  (End  Plates  Off), 
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Figure  68.  Flat  Plate  Oil  Flow  Study  (End  Plates  Off), 


Figure  69.  Flat  Plate  Oil  Flow  Study  (End  Plates  Off), 


Figure  70.  Flat  Plate  Oil  Flow  Study  (End  Plates 


Figure  71 .  Flat  Plate  Oil  Flow  Study  (End  Plates 
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Figure  72.  Effect  of  End  Plates  on  Pressure  Distribution,  6  =  10°.  M  s  k 
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Figure  73.  Effect  of  End  Plates  on  Pressure  Distribution,  fi  =  20°,  M 
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No  extensive  transverse  effects  of  the  endplates  on  pressure  are  noted,  with 
the  exception  of  the  case  with  a  =  0°,  8  =  10°  and  Rero/ft  *  1.1  x  10^.  The 
characteristic  pressure  distribution  interact: on  distances  in  the  longitudi¬ 
nal  direction  remained  approximately  the  same  for  end  plates  "on"  or  "off" . 

The  transverse  pressure  distribution  in  the  re-attachment  region  is  strongly 
influenced  by  decreasing  Reynolds  number  and  shows  local  peaks  at  both  the 
center  of  the  plate  and  toward  the  sides ,  at  least  as  far  as  the  instrumenta¬ 
tion  extends  (Figures  72  and  73).  This  trend  is  particularly  well  shown  in 
Figure  72  at  Re^/ft  =  1.1  x  10°.  It  appears  that  a  thick  boundary  layer  at 
low  Reynolds  numbers  promotes  the  three-dimensionality  of  the  flow.  In 
general,  however,  the  transverse  pressure  distribution  upstream  of  the  hinge 
lino  is  farily  insensitive  to  Reynolds  number  (Figures  72  and  73). 

6 

For  Mo,  =  8  and  Re^/ft  =  3.3  x  10  tests  were  made  at  a  =  0°  and  15°  for 
5f=  10c  ,  20° ,  and  30°  (Reference  59) •  It  is  observed  that  at  the  given  flow 
conditions  for  a  =  0°  the  pressure  reached  on  the  flap  is  higher  for  the 
end  plates  off  tests  (Figures  7k  and  75).  Ato=  15°  and  6f  =  30°  the  flap 
pressures  seem  independent  of  end  plate  pressure  (Figure  75).  This  is  con¬ 
sistent  vrith  the  observations  at  Re^/ft  =  6.6  x  10°,  as  the  compression  angle 

of  attack  increased  the  Ipcal  Reynolds  number  ahead  of  the  flap  from 
3.3  x  10°/ft  to  5.88  x  10°/ ft.  Extensive  pressure  variation  in  the  spanvise 
direction  is  found  both  upstream  and  downstream  of  the  hinge  line.  For  a  =  0°  , 
the  upstream  pressure  distribution  is  generally  consistent  with  the  =  5  test 
results .  The  plateau  region  exhibits  slightly  higher  pressures  with  the  end 
plates  on  (Figures  7k  and  75).  No  significant  change  in  the  upstream  and 
downstream  parameters  is  observed. 

DELTA  WING  CONFIGURATIONS 


The  generalized  correlations  for  shockwave-boundary  layer  interaction  were 
derived  from  experimental  data  obtained  for  models  with  essentially  two- 
dimensional  flow:  flat  plate-ramp,  flat  plate  with  incident  shock,  and  flat 
plate  vrith  a  step.  To  examine  the  general  applicability  of  these  results, 
the  flow  over  various  delta  configurations  is  discussed. 

Flow  About  a  Delta  Wing 


The  flow  about  a  delta  wing  configuration  is  rather  complex;  however,  a  number 
of  investigations  on  delta  shapes  has  been  conducted  (References  62  through  66) . 
A  very  useful  method  of  studying  comp], ex  flow  is  the  oil  now  method .  Figure 
76  shows  a  surface  flow  pattern  inferred  from  oil  flow  photographs  on  the  com¬ 
pression  surface  of  a  sharp  leading  edge  delta  wing  at  various  angles  of  in¬ 
cidence  (Reference  66) .  The  surface  flew  and  the  boundary  layer  flow  at  small 
angles  of  attack  have  an  inward  component  which  is  attributed  to  the  boundary 
layer  induced  pressure  gradient  at  the  leading  edge.  Conversely,  the  external 
flow  is  directed  away  from  the  centerline.  At  higher  angles  of  attack  the 
flow  turns  parallel  and  eventually  diverges .  At  an  incidence  of  about  10°  the 
flow  divides  along  a  line  at  about  90$  of  the  span,  with  the  outooard  flow- 
directed  toward  the  leading  edge.  At  higher  angle  of  attack,  the  dividing  line 
moves  inboard  and  reaches  the  centerline  of  the  body  at  a  s-  18°  Beyond 
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Figure  75.  Effect  of  End  Plates  on  Pressure  Distribution,  6^=  30°,  Mn=  8 
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a  s  20"  the  flow  field  is  such  that  the  leading  edge  of  the  delta  wing  may 
become  a  trailing  edge.  This  outboard  flow  is  a  result  of  shock  losses  which 
are  predominant  at  higher  incidence. 

Plow  patterns  on  the  lee  surfaces  of  delta  wings  at  various  angles  of  attack 
and  Mach  numbers  have  been  examined  in  Reference  62.  At  small  angles  of  in¬ 
cidence  the  flow  was  noted  to  be  attached  over  the  entire  upper  surface.  This 
is  observed  in  Figure  77  for  lower  Mach  numbers  and  in  Figure  77  at  higher 
Mach  numbers. 

As  the  angle  of  attack  increases  for  the  low  Mach  number  case  the  flow  separ¬ 
ates  from  the  leading  edge.  A  coiled  vortex  sheet  is  thus  formed  which  lies 
across  the  expansion  surface  along  the  leading  edges  (Figure  77).  A  span- 
wise  outflow  is  induced  beneath  the  vortex  and  a  secondary  separation  occurs 
along  the  leading  edge  (Figure  77).  At  lower  sweep  angles  the  vortex  sheet 
re-attaches,  forming  a  long  separation  bubble. 

For  higher  Mach  numbers,  the  flow  is  attached  at  lower  angles  of  attack.  A 
weak  shock  is  noted  to  emanate  from  the  vortex  region  and  deflect  the  flow. 

When  this  shock  attains  sufficient  strength  the  boundary  layer  separates  and 
forms  a  coiled  vortex  sheet  (Figure  77).  Again,  at  lower  angles  of  attack, 
the  boundary  layer  re-attaches  to  the  surface  and  forms  a  separation  bubble. 

On  the  expansion  surface  of  a  thick  delta  wing  at  incidences  at  which  suction 
is  developed,  separation  lines  become  noticeable  and  moves  outboard  as  the 
incidence  is  increased  (Figure  78).  The  area  between  the  separation  lines  on 
both  sides  of  the  centerline  is  an  area  of  reduced  shear  and  thus  the  flow  is 
inclined  awray  from  the  centerline.  At  a  >  25°  re-attachment  of  the  separated 
flow  occurs  at  the  centerline  of  the  body. 

At  a  s  19°  there  is  visible,  a  concentration  of  oil  filaments  from  the  turning 
of  the  flow  away  from  the  plane  of  symmetry.  This  indicates  a  flow  retardation 
due  to  a  presence  of  .an  external  shock  wave.  Deviation  from  conical  flow  is 
seen  from  the  curvature  of  the  separation  line  near  the  trailing  edge.  This 
is  attributed  to  the  upstream  propagation  through  the  thick  laminar  boundary 
layer  of  strong  pressure  gradients  generated  by  compression  at  the  trailing 
edge. 

The  flow  structure  of  the  resulting  conical  flow  is  chown  in.  the  spanwise  cross- 
section  in  Figure  79*  This  figure  identifies  the  location  of  the  bow  shock 
and  expansion  flow  field  emanating  from  the  leading  edge.  The  Doundary  layer 
is  thick  and  the  adverse  pressure  gradient  due  to  external  compression  is 
propagated  toward  the  leading  edge . 

The  position  of  separation  determined  from  the  oilflow  photographs  is  shown  in 
Figure  80  (Reference  66).  Although  flow  details  near  the  centerline  are  lack¬ 
ing,  it  is  suggested  that  rolling  up  of  the  separated  blow  promotes  a  new 
expansion  of  the  external  flow  which  then  turns  parallel  to  the  surface  and 
gives  rise  to  compression  at  the  centerline.  An  idealized  flow  over  such 
configurations  has  been  analyzed  and  verified  by  Fowell  (Reference  63) .  It 
was  shown  that  a  Prandtl -Meyer  expansion  centered  at  the  leading  edge  turns 
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the  flow  inward  and  at  some  distance  inboard  recompression  takes  place  through 
an  oblique  shock  wave  which  turns  the  flow  parallel.  Considering  viscous 
effects  in  a  real  flow,  a  strong  recompression  may  cause  flow  separation  and 
formation  of  a  vortex.  The  flow  outside  o'f  the  vortex  may  reattach,  and 
inboard  of  the  reattachment,  flow  is  turned  parallel  to  the  centerline. 

The  flow  patterns  exhibited  by  the  blunt  delta  wing  are  roughly  similar  to 
those  presented  for  the  sharp  leading  edge  (Reference  65) .  The  flow  pattern, 
at  a  *  0°  is  characterized  by  central  outflow  as  a  result  of  high  induced 

pressure  at  the  apex,  and  is  contained  by  the  shock  and  the  induced  high  pres¬ 
sure  along  the  leading  edges  (Figure  81) .  Surface  flow  on  cylindrical  edges 
is  inward  and  follows  the  pressure  gradient.  At  increasing  angle  of  attack  an 
inward  flow  develops  until  a  >_  20°  where  an  outward  flow  occurs  (Figure  81) . 

An  investigation  of  blunted  delta  wings  with  deflected  flaps  is  presented  in 
Reference  67.  The  existence  of  deflected  flaps  on  the  model  surface  modifies 
the  flow  over  a  substantial  portion  of  the  configuration.  The  oil  film 
patterns  for  this  model  disclose  a  very  complex  flow  field.  The  surface  flow 
on  the  compression  surface  shows  a  strong  influence  of  the  blunt  apex  in  the 
central  area  of  the  model  (Figure  82) .  A  definite  "blast"  region  clear  of 
oil  is  observed  in  the  central  high  pressure  area.  Separation  lines  upstream 
of  both  flaps  are  well  developed  and  symmetrical.  A  flow  pattern  at  a  *  0° 

is  shown  in  Figure  83.  Attached  flow  exists  near  the  leading  edges.  Further 
inboard,  flow  separation  occurs  with  reattachment  along  the  centerline.  Fiow 
separation  upstream  of  the  flaps  is  well  outlined.  Separation  effects  pro¬ 
pagate  upstream  and  form  complex  flow  patterns  a  significant  distance  ahead  of 
the  flaps.  Similar  flow  with  more  details  exists  on  the  expansion  surface  at 
a  =  -30°  (Figure  84) . 

Pressure  Distribution  on  Flat  Delta  Wing 

The  general  applicability  of  the  two-dimensional  correlations  for  separated 
flow  to  more  complex  configurations  will  now  be  investigated. 

A  blunted  flat-bottom  delta  wing  from  Reference  67  is  considered  first  (Figure 
85).  In  order  to  verify  the  existence  of  two-dimensional  flo*',  the  experimental 
pressure  distribution  in  the  streamwise  and  spanwise  directions  is  considered. 
For  the  flat-bottom  delta  configuration,  the  streamwise  and  spanwise  experi¬ 
mental  distribution  at  33°  angle,  of  attack  (M»  *  8,  A  *  70°)  is  shown  in 

Figures  86  and  87.  It  is  noted  that  the  pressure  distribution  downstream  of 
the  nose  agrees  very  closely  with  the  two-dimensional  oblique  shock  value 
obtained  from  Reference  36.  This  indicates  that  conditions  in  the  separated 
region  upstream  of  the  flaps  may  be  determined  by  using  two-dimensional 
correlations. 

The  local  Mach  number  in  the  vicinity  of  the  flap  was  determined  from  the  ratio 
of  flap  pressure  P£-^a_  to  undisturbed  pressure,  Pa.  This  was  accomplished 
by  finding  M^in  0,  where  0  is  the  shock  angle,  through  the  use  of 
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Figure  82.  Oil  Film  Photograph  of  Lower  Surface  of  Blunt  Pyramidal  Configuration 
(Ref.  67);  Mc0=  5,  «  =  14.3°  6f  =  40° 
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Figure  84.  Oil  Film  Photograph  of  Lower  Surface  of  Blunt  Pyramidal  Configuration 
(Ref .  67);  Mro  =  5,  a  =  -30°,  6f  =  40 
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X'  =  NONDIMEN  SIGNAL  STREAMWISE  DISTANCE  FROM  VIRTUAL  APEX 


Figure  86.  Lengthwise  Pressure  Distribution  Over  Pyramidia'l  Configuration 
(Flat  Bottom  Surface) 
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OBLIQUE  SHOCK  METHOD 


-/-.a£  in  the  normal  shock  tables  of  Reference  36.  is  then  obtained  by 
P  W 

iteration  with  the  aid  of  Chart  II  of  the  same  reference .  The  value  for  M  a 

was  found  to  be  almost  identical  to  the  two-dimensional  value  corresponding 
to  the  given  angle  of  attack.  Therefore  the  local  Reynolds  number  is  found 
through  the  use  of  the  oblique  shock  tables  of  Reference  36. 


The  distance  and  pressure  parameters  d-^ ,  if^,  Cpp?  aiI<i  Cpg  flow 

upstream  of  the  flaps  were  calculated  as  noted  in  Section  IV.  The  existence 
of  laminar  flow  (Re  =  7.81  x  105.  Mq.  =  4.79)  was  verified  by  the  Deem 

T>  . 

Murphy  formula  v Equation  17) . 


Reasonable  agreement  with  the  experimental  values  of  Reference  67  was  obtained 
(Figures  88  through  91).  Since  the  Re^changes  in  the  spanwise  direction, 


the  calculated  parameters 
mended  that  an  average  Re 


will  be  a  function  of  y. 
^  be  taken  to  determine 
^0 


It  is  therefore  recom- 
the  pressure  profile.  For 


small  values  of  A  and  for  downstream  flap  positions  this  will  not  be  necessary. 
In  cases  where  the  control  surfaces  are  located  in  the  region  of  influence  of 
the  blunt  leading  edge,  bl  witness  effects  must  be  considered  for  moderate  angle 
of  attack.  The  blast  wave  analysis  of  Reference  34  may  be  used  to  find  the 
pressure  increment  due  to  the  blunt  leading  edge.  Therefore, 
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The  mviscid  pressure  distribution  is  thus: 
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This  may  be  solved  for  Pq/Pq,  at  each  xQ  position  .  If  viscous  effects  are 
found  to  be  important  by  the  previously  mentioned  criteria  (Section  III), 
they  may  be  accounted  for  by  the  methods  presented  in  that  section.  The 
preceding  analysis  assumes  the  flow  to  be  two  dimensional  in  the  sense  that 
strip  theory  is  applicable. 

Flow  About  a  Delta  Wing  with  Dihedral 

The  dihedral  surface  of  the  pyramidal  configuration  (Figure  83)  exhibits  a 
significantly  different  flow  pattern  from  that  of  the  flat  bottom  surface. 
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Figure  88.  Pressure  Distribution  Over  Pyramidal  Configuration  (Fiar  Bottom  Surface) 
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Figure  89.  Pressure  Distribution  Over  Pyramidal  Configuration 
(Flat  Bottom  Surface),  Re^3  3.3  x  loVft 
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Figure  9G.  Pressure  Distribution  Over  Pyramidal  Configuration  (Flat  Bottom 
Surface),  X/. *  =  0.75,  R  =  1  .1  X  10  /ft,  8,  =  *0° 
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Figure  91 .  Pressure  Distribution  Over  Pyramidal  Configuration  (Flat  Bottom  Surface) 
8f  =  «• ,  X^,  =  0.85.Reco  =  1 .1  x  106/ft 
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The  oil  flow  pattern  (Figures  4  through  7,  Reference  67)  indicates  outboard 
flow  on  both  sides  of  the  ridge.  A  well  defined  blast  region  no  longer 
exists .  The  interaction  region  ahead  of  the  flap  undergoes  a  severe  dis¬ 
tortion  in  the  transverse  direction .  The  separated  region  is  swept  outward 
and  is  noted  to  flow  over  the  blunt  edges. 

An  outflow  is  seen  to  exist  all  along  the  leading  edges  for  all  angles  of 
attack  at  which  the  pressures  on  the  dihedral  surfaces  exceed  those  on  the 
lower  surface .  Consequently,  the  resultant  flow  field  around  the  model  is 
rather  complex  and  the  application  of  certain  approximations  in  calculating 
the  flow  properties  seems  well  justified.  The  accuracy  of  these  approxima¬ 
tions  will  be  established  by  comparison  with  experimental  values . 

Pressure'  Distribution  on  a  Delta  Wing  with  Dihedral 


To  test  the  applicability  of  the  two-dimensional  analysis  to  a  blunt  delta 
wing  (A=  70°)  with  dihedral,  the  local  flow  conditions  on  dihedral  surfaces 
were  examined.  It  was  found  that  the  pressure  in  the  streamwise  direction  is 
almost  constant  and  decreases  only  slightly  in  the  spanwise  direction  as  shown 
on  Figures  $2  and  93  for  a  =  29°  ;  conical  flow  may  be  more  applicable  as  con 
eluded  from  the  conical  shock  shape  appearing  at  higher  angles  of  attack 
(Reference  68) .  The  instrumentation  density  in  the  interaction  region  was  not 
sufficient  to  allow  the  determination  of  the  characteristic  parameters  dl>  if* 
and  Cpp. 
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(Dihedral  Surface) 


Section  VI 


PREDICTION  OF  EFFECTIVENESS  OF  AERODYNAMIC  CONTROLS 


This  section  develops  a  systematic  procedure  for  the  application  of  correla¬ 
tion  results  presented  in  preceding  sections  to  the  prediction  of  hypersonic 
control  effectiveness.  Wherever  practicable,  design  charts  are  developed  to 
facilitate  the  calculation  procedure.  The  methods  are  intended  for  use  over 
the  broad  range  of  flight  conditions  encountered  by  a  typical  hypersonic 
cruise/glide  vehicle.  It  is  to  be  enphasized,  however,  that  the  methods  are 
based  upon  data  obtained  for  a  finite  range  of  experimental  conditions  and 
that  the  validity  of  an  extrapolation  of  the  methods  for  conditions  far  out¬ 
side  the  experimental  range  is  yet  to  be  established.  The  range  of  variables 
covered  by  the  experimental  data  is  indicated  in  the  correlation  figures 
(Section  IV)  and  prediction  charts.  Application  of  the  procedure  described 
in  this  section  for  calculation  of  control  effectiveness  is  illustrated  in 
the  sample  problems  presented  in  Appendix  II. 

FLIGHT  CONDITIONS 


The  altitude-velocity  envelope  of  a  vehicle  defines  the  range  of  free-stream 
conditions  which  are  encountered  during  flight.  For  a  given  velocity,  the 
maximum  altitude  at  which  a  vehicle  can  fly  is  determined  by  the  minimum 
dynamic  pressure  which  will  sustain  flight,  and  the  minimum  altitude  is  re¬ 
stricted  by  the  capability  of  a  vehicle  to  withstand  aerodynamic  heating. 

The  aerodynamic  and  heating  boundary  curves  for  a  typical  hypersonic  cruise/ 
glide  vehicle  were  obtained  from  Reference  50  and  are  presented  in  Figure  9^ 
The  information  developed  in  this  section  is  confined  to  freestream  flight 
conditions  in  the  altitude-velocity  envelope  for  Mach  numbers  greater  than  5* 

The  angle  of  attack  at  which  a  vehicle  flies  is  limited  by  vehicle  performance 
and  structural  heating  considerations .  The  angle  of  attack  below  which  the 
lift/drag  ratio  5s  favorable  and  structural  heating  is  not  excessive  is  ap¬ 
proximately  20  degrees.  Since  this  study  is  concerned  primarily  with  flat 
plate  surfaces,  angle  of  attack  also  represents  flow  deflection  angle. 
Therefore,  angle  of  attack  is  used  to  denote  flow  deflection  angle  and  is 
limited  to  the  range  from  0  to  20  degrees. 

The  range  of  control  surface  deflection  angle  from  0  to  30  degrees  relative 
to  the  forward  surface  is  assumed  to  encompass  most  cases  of  practical 
interest. 
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Figure  94.  Typical  Flight  Envelope  for  Lifting  Hypersonic  Vehicles 
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LOCAL  FLOW  CONDITIONS 


Local  flow  properties  are  u^ed  to  define  the  parameters  .-/hich  describe  com¬ 
pression  corner  flow  separation.  The  first  problem,  then,  in  a  control  effec¬ 
tiveness  analysis  is  the  determination  of  flow  conditions  upstream  of  the 
control  surface.  Section  III  provides  a  discussion  of  methods  for  predicting 
local  flow  properties  along  blunted,  planar  surfaces.  It  was  shown  that  on 
windward  surfaces  far  downstream  of  bhe  leading  edge,  local  flow  properties 
are  approximately  the  same  as  flow  conditions  behind  an  oblique  shock. 

For  many  practical  reentry  configurations,  control  surfaces  are  located  near 
the  trailing  edge  of  essentially  planar  surfaces  and  are  sufficiently  far 
from  the  leading  edge  for  bluntness  effects  to  be  negligible.  Also,  since 
leeward  surfaces  in  hypersonic  flow  are  generally  ineffective,  the  problem  is 
reduced  to  investigation  of  compression  corner  flow  separation  on  windward 
surfaces.  Therefore,  in  many  cases  the  flow  properties  forward  of  a  control 
surface  can  be  assumed  to  be  approximately  the  same  as  conditions  behind  an 
oblique  shock. 

Since  oblique  shock  assumptions  are  often  valid,  oblique  shock  properties 
have  been  determined  for  a  number  of  flight  conditions  in  the  altitude- 
velocity  envelope  (Figure  9W .  Section  3-4  recommends  that  either  perfect 
gas  or  equilibrium  real  gas  assumptions  be  used  for  an  oblique  shock  depend¬ 
ing  on  velocity-altitude  conditions  and  flow  deflection  angle. 

For  the  flight  regimes  where  fJow  properties  deviate  from  those  given  by 
oblique  shock  theory  for  a  perfect  gas,  equilibrium  real  gas  effects  were 
determined  by  means  of  similarity  parameters  from  approximate  theory  and 
effective  ratio  of  specific  heats  presented  in  Reference  69.  This  reference 
presents  the  effective  ratio  of  specific  heats,  7  ,  at  various  altitudes 
(based  on  ARDC  1959  model  atmosphere)  as  a  function  of  the  component  of  the 
Mach  number  normal  to  the  shock  wave.  The  solution  for  the  correct  7e  to  be 
used  in  the  calculations  involved  an  iterative  procedure  for  determining  the 
shock  wave  angle  and  corresponding  normal  component  of  the  Mach  number.  For 
the  first  iteration,  the  Mach  number  component  normal  to  a  perfect  gas  ob¬ 
lique  shock  angle  was  assumed  for  determining  the  initial  7  (Reference  69, 
Figure  12).  This  value  of  7p  was  used  to  calculate  the  value  of  the  param¬ 
eter  (7  +  l)  sin  a.  Figure  4  of  Reference  69  presents  curves  which  cor- 
e  p 

relate  the  normal  Mach  number  component  as  a  function  of  this  parameter.  The 
normal  Mach  number  obtained  from  these  curves  was  used  for  the  second  itera¬ 
tion.  This  procedure  was  repeated  until  the  initial  and  final  normal  Mach 
number  components  converged. 

Having  determined  the  effective  ratio  of  specific  heats,  the  oblique  shock 
correlations  of  Reference  69  were  used  for  rapid  calculation  of  the  flow 
properties  for  an  equilibrium  real  gas. 
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The  oblique  shock  Reynolds  number  ratio  was  determined  from  the  following 
expression: 


P  u 

O'  Q1 

P»u« 


) 


O.67 


This  expression  is  derived  using  an  exponential  viscosity-temperature  relation 
with  a  power  of  0.6?«  Comparison  of  the  exponential  vi scosity- temperature  re¬ 
lation  for  various  powers  with  Sutherland's  formula  is  shown  in  Figure  95* 
This  comparison  indicates  that  over  the  range  of  temperature  of  practical  in¬ 
terest,  the  O.67  power  provides  sufficient  accuracy  for  this  application. 
Although  a  lower  power  provides  better  agreement  at  high  temperature,  the 
disagreement  is  increased  in  the  intermediate  range  of  temperature  where  most 
of  the  design  conditions  are  expected  to  occur. 


The  Mach  number  behind  an  oblique  shock  for  a  real  gas  was  evaluated  with  the 
aid  of  References  69  and  70  using  the  relation 
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The  oblique  shock  velocity  ratio  was  determined  from  equation  25a  of  Reference  69. 
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The  speed  of  sound  ratio  was  obtained  from  data  presented  in  Table  IV  of 
Reference  70. 


Figures  96  through  99  present  the  pressure,  temperature,  Mach  number,  and 
Reynolds  number  behind  an  oblique  shock  for  flow  deflection  angles  from  0  to 
20  degrees  for  the  following  altitude  and  Mach  number  conditions: 


Altitude  (feet) 

50,000 

100,000 

150,000 

200,000 

250,000 

300,000 


Mach  Number 

5,  7 
5,  7,  10 
10,  15 
15,  20,  30 
20,  30 
20,  30 


The  selection  of  Mach  number  is  restricted  to  those  Mach  numbers  in  the 
vicinity  of  the  flight  envelope  for  which  Reference  69  presents  data. 
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Figure  95.  Comparison  of  Sutherland's  Viscosity  Law  with  Exponential  Viscosity  - 
Temperature  Relation 


ANGLE  OF  INCIDENCE,  a  ,  DEGREES 
(a)  ALTITUDE  =50,000  FT 


Figure  96.  Oblique  Shock  Pressure  Ratio  Including  Real  Gas  Effects 
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(b)  ALTITUDE  =  100,000  FT. 


Figure  96.  Oblique  Shock  Pressure  Ratio  Including  Rea!  Gas  Effects  (Continued) 
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ANGLE  OF  INCIDENCE,  a,  DEGREES 
(c)  ALTITUDE  =  150,000  FT. 


Figure  96.  Oblique  Shock  Pressure  Ratio  Including  Real  Gas  Effects  (Continued) 


o  5  10  15  20 

ANGLE  OF  INCIDENCE,  a  ,  DEGREES 
(f)  ALTITUDE  =  300,000  FT, 

Figure  96.  Oblique  Shock  Pressure  Ratio  Including  Real  Gas  Effects  (Concluded) 


(a)  ALTITUDE  =  50,000  FT. 


Figure  97.  Oblique  Shock  Temperature  Ratio  Including  Real  Gas  Effects 


158 


0  5  10  15  20 

ANGLE  OF  INCIDENCE, a  ,  DEGREES 


(c)  ALTITUDE  =  150,000  FT. 


Figure  97.  Oblique  Shock  Temperature  Ratio  Including  Real  Gas  Effects  (Continued) 
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ANGLE  OF  INCIDENCE,  a  ,  DEGREES 
(d)  ALTITUDE  =  200,000  FT. 


Figure  97.  Oblique  Shock  Temperature  Ratio  Including  Real  Gas  Effects  (Continued) 


l6l 


PERFECT  GAS  | 

EQUILIBRIUM  REAL  GAS 


0  5  10  15  20 

ANGLE  Or  INCIDENCE,  a  ,  DEGREES 

(e)  ALTITUDE  =250,000  FT. 
igure  97.  Oblique  Shock  Temperature  Ratio  Including  Real  Gas  Effects  (Continued) 
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ANGLE  OF  INCIDENCE,  a  ,  DEGREES 
(f)  ALTITUDE  =  300,000  FT. 

Figure  97.  Oblique  Shock  Temperature  Ratio  Including  Real  Gas  Effects  (Concluded) 
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Figure 98.  Mach  Number  Behind  an  Oblique  Shock  Including  Real  Gas  Effects 
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ANGLE  OF  INCIDENCE,  a  ,  DEGREES 
(b)  ALTITUDE  =  100,000  FT. 

Figure  98.  Mach  Number  Behind  an  Oblique  Shock  Including  Real  Gas  Effects  (Continued) 
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ANGLE  OF  INCIDENCE,  a  ,  DEGREES 
(c)  ALTITUDE  =  150,000  FT. 

Figure  98.  Mach  Number  behind  an  Oblique  Shock  Including  Real  Gas  Effects  (Continued) 


166 


167 


ANGLE  OF  INCIDENCE,®  ,  DEGREES 
(a)  ALTITUDE  =50,000  FT. 


Figure  99.  Oblique  Shock  Reynolds  Number  Ratio  Including  Real  Gas  Effects 
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(b)  ALTITUDE  =  100,000  FT. 


Figure  99.  Oblique  Shock  Reynolds  Number  Ratio  Including  Real  Gas  Effects  (Continued) 


ANGLE  OF  INCIDENCE,  «  ,  DEGREES 
(c)  ALTITUDE  =  150,000 


Figure  99.  Oblique  Shock  Reynolds  Number  Ratio  Including  Real  Gas  Effects  (Continued) 
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ANGLE  OF  INCIDENCE,  «  ,  DEGREES 
(d)  ALTITUDE  =  200,000  FT. 


Figure  99. 


Oblique  Shock  Reynolds  Number  Rario  Including  Real  Gas  Effects  (Continued) 
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Figure  99. 


Oblique  Shock  Reynolds  Number  RaHo  Including  Real  Gas  Effects  (Continued) 
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Figure  99.  Oblique  Shock  Reynolds  Number  Ratio  Including  Real  Gc;  Effects  (Conductor! 


BOUNDARY  LAYER  THICKNESS 


Since  boundary  layer  thickness  is  the  parameter  used  for  non-dimens ionali zing 
the  interaction  distances  of  a  separated  flow  field,  charts  have  been  con¬ 
structed  to  simplify  the  procedure  for  calculating  boundary  layer  thickness. 

Section  III  shows  that  hypersonic  boundary  layer  thickness  may  be  determined 
using  Reynolds  number  obtained  by  the  reference  temperature  method  (Ref.  71) . 
Using  the  0.67  power  viscosity-temperature  relation,  the  ratio  of  reference 
temperature  Reynolds  number  to  local  Reynolds  number  reduces  to  the  following 
form: 


Substitution  of  the  equations  which  define  reference  temperature  and  adiabatic 
wall  temperature  into  the  above  equation  results  in  the  expression 


Re* 


a 


T 

0.28  +  0.5  ~  +  0.22 

a 


1  + 


M  2Prn 
or 


-1.67 


Reference  temperature  Reynolds  number  calculations  were  based  upon  a  Prandti 
number  of  0.72  with  n  equal  to  l/2  for  laminar  flow  and  to  l/3  for  turbulent 
flow.  Wall  temperatures  were  restricted  to  the  expected  range  of  practical 
interest  (l<Tw/Ta  <  IO^T^/Tq.  ).  Figures  100  and  101  present  the  reference 
temperature  Reynolds  number  ratio  for  laminar  and  turbulent  flow  as  a  function 
of  local  Mach  number  for  various  ratios  of  wall  to  local  stream  temperature. 

Plots  of  the  laminar  and  turbulent  boundary  layer  thickness  equations,  based 
on  reference  temperature  Reynolds  number,  are  presented  i*.  Figures  102  and  103. 
Thr  equations  shown  omit  the  pressure  correction  factor  which  should  be  in¬ 
cluded  in  the  case  of  significant  viscous  interaction  or  bluntness  induced 
pressure  gradient  (See  Section  III). 

PRESSURE  DISTRIBUTION 


The  development  of  semi-empirical  correlations  for  describing  the  pressure 
distribution  in  the  laminar  and  turbulent  separated  region  was  discussed  in 
Section  IV.  The  pertinent  expressions  have  been  presented  in  the  form  of  work¬ 
ing  charts  to  facilitate  calculation  of  the  surface  pressure  distribution  in 
the  region  of  separated  compression  corner  flow.  The  charts  are  shown  in 
Figures  104  through  114.  The  procedure  for  using  the  charts  is  discussed  be¬ 
low,  ard  their  use  is  illustrated  in  the  sample  calculations  presented  in 
Appendix  II. 
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Figure  101 .  The  Effects  of  Mach  Number  and  Wall  Temperature  on  Local  Reference 
Temoerature  Reynolds  Number  Ratio  (Turbulent) 
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Figure  1 02 .  Laminar  Boundary  Layer  Thickness 
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Figure  104.  Pressure  Coefficient  for  Incipient  Boundary  Layer  Separation  in  Laminar  Flow 
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Figure  105.  Pressure  Coefficient  for  Incipient  Boundary  Layer  Separation  in  Turbulent  Flow 


Figure  106.  Laminar  Plateau  Preaure  Coefficient 
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Figure  10 7,  Turbulent  Plateau  Pressure  Coefficient 
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Figure  108.  Upstream  Interaction  Distance  for  Laminar  Flow 
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Figure  110.  Free  Interaction  Length  for  Laminar  Flow 


Figure  112.  Downstream  Interaction  Distance  to  Peak  Pressure  on  Flap  for  Laminar  Flow,  d; 


LIMITING  CONDITIONS: 


1.  IF-pSl,  THEN 
dl 


d. 


!.  IF ^-<0.25,  THEN~^==  =2  — 

d 


1 


(Cf  =  FLAP  CHORD) 


Figure  113.  Downstream  Interaction  Distance  to  Peak  Pressure  for  on  Flap  Turbulent  Flow,  d ^ 
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Figure  1 14.  Downstream  interaction  Distance  to  Pressure  Rise  on  Flap  for  Laminar  Flow, 
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The  first  step  in  defining  the  pressure  distribution  is  to  determine  if  flow 
separation  exists.  The  pressure  coefficient  data  recommended  for  predicting 
incipient  boundary  layer  separation  in  laminar  or  turbulent  flow  are  presented 
in  Figure  100. 

It  has  been  verified  (Section  4.1)  that  the  existence  of  laminar  or  turbulent 
flov  can  be  predicted  with  adequate  accuracy  by  the  Deem-Murphy  criterion 
(Reference  51  )> 

If  flow  separation  exists,  an  iterative  procedure  is  required  to  determine 
the  separation  location.  For  a  given  local  Mach  number,  the  upstream  inter¬ 
action  distance  is  a  function  of  the  plateau  pressure  level,  which  is  a  func¬ 
tion  of  the  Reynolds  number  at  the  separation  point.  (The  interaction  flow 
model  is  discussed  in  Section  IV  and  is  shown  in  Figure  28.)  Therefore, 
for  any  separation  location  a  corresponding  plateau  pressure  level  and  up¬ 
stream  interaction  distance  is  defined.  In  solving  for  the  separation  loca¬ 
tion,  the  upstream  interaction  distances  corresponding  to  a  number  of  assumed 
separation  locations  are  calculated  until  the  point'  is  found  for  which  the 
sum  of  the  downstream  distance  to  the  separation  point  and  the  upstream  in¬ 
teraction  distance  is  equal  to  the  distance  to  the  hinge  line.  The  procedure 
is  as  follows:  * 

1.  Assume  a  separation  location  Xq  (the  prime  superscript  denotes 
distance  measured  from  an  effective  starting  point  for  equiva¬ 
lent  flat  plate  boundary  layer  generation),  and  calculate  the 
Reynolds  number. 

2.  Read  the  plateau  pressure  coefficient  from  Figure  106  or  107. 

3.  Obtain  the  upstream  interaction  distance,  d^,  non-dimensionalized 
by  boundary  layer  thickness,  from  Figure  108  or  109* 

4.  For  a  given  wall  temperature,  calculate  the  reference  temperature 
Reynolds  number  from  Figure  100  or  101  and  the  corresponding  bound¬ 
ary-layer  thickness  from  Figure  102  or  103. 

5.  Compute  d^  and  Xq  +  d^.  If  Xq  +  d^  ^  X^,  repeat  procedure. 

Figures  no  through  114  present  the  remaining  interaction  distances  (£fi>  d2> 
and  d^)  which  define  the  pressure  distribution,  as  functions  of  the  plateau 
pressure  level,  upstream  interaction  distance,  and  boundary  layer  thickness. 


INCREMENTAL  FORCE  AND MOMENT  EQUATIONS 


Equations  have  been  derived  for  the  two-dimensional  force  and  moment  incre¬ 
ments  produced  by  separated  compression  corner  flow.  The  force  and  moment 
equations  were  developed  in  terms  of  the  distances  and  pressure  levels  which 
define  the  pressure  distribution.  A  body  fixed  coordinate  axes  system  with  ori¬ 
gin  at  an  arbitrary  vehicle  moment  reference  point  (MRP)  was  used  in  the 
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derivation.  The  coordinate  axes  system  and  convention  foi  positive  sense  of  the 
forces  and  moments  are  shown  in  Figure  115.  The  equations  are  summarized  below. 
The  complete  derivation  of  the  equations  is  presented  in  Appendix  I.  It  is  to 
be  emphasized  that  the  following  equations  represent  the  force  and  moment  incre¬ 
ments  produced  on  the  windward  surface  of  a  deflected  control  and  that  incre¬ 
ments  produced  on  the  leeward  surface  are  not  included;  the  subscript  w  is  used 
to  denote  this  fact.  In  many  cases,  however,  the  leeward  contribution  is  negli¬ 
gible.  Also,  the  effect  of  skin  friction  control  effectiveness  was  not  included. 
Note  that  all  force  increments  and  pitching  moment  increments  are  based  on  the 
pressure  on  the  windward  side  at  the  beginning  of  interaction.  Windward  compo¬ 
nent  of  elevon  hinge  moment  is  based  on  absolute  pressure  distribution. 

1.  Normal  force  increment: 


2.  Increment  of  pitching  moment  about  HL  due  to  normal  force  increment: 


3-  Center  of  pressure  of  normal  force  increment: 


AC  L 

ro-HL 


AN 


All. 


w 


AC. 


w 


N. 


w 


4,  Axial  force  increment: 


(38) 


M  , 

t> 


+ 


(39) 
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figure  115.  Two-Dimensional  Force  and  Moment  Convention  for  Symmetrical  Elevon  Deflection 


5.  Increment  of  pitching  moment  about  HL  due  to  axial  force  increment: 


<CP  J, 


(c  ) 

'  poa 


Jl 


2  2 

C  -  d  ^  -  d^d 
f  2 

2  T" 


F'-f 


]} 


6.  Center  of  pressure  of  axial  force  increment: 


Z 

cPaa 

v 


(40) 


(41) 


7.  Increment  of  pitching  moment  about  vehicle  MRP: 


%y. 


=  AC„.  X  •  -  AC.  Z 

N  crA*T  A,  cp  .  . 

w  *AN  w  ^  AA. 

w  w 


(42) 


8.  Windward  component  of  elevon  hinge  moment  (based  on  absolute  pressure) 


HM 
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-  (c  )  c  2  -  1.  [(< 
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2  2 
d..  +  d„d„  +  d  ^ 
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These  equations  provide  a  direct  means  of  calculating  the  two-dimensional, 
windward  effectiveness  of  a  control  surface  once  the  idealized  pressure  dis¬ 
tribution  for  compression  corner  flow  is  known. 

COMPARISON  OF  EXPERIMENTS  WITH  PREDICTIONS 

To  illustrate  the  applicability  of  the  prediction  methods  for  estimating  con¬ 
trol  effectiveness,  incremental  force  and  moment  coefficients  have  been  cal¬ 
culated  for  a  test  configuration  and  compared  with  the  experimental  data. 
Reference  72  presents  force  test  results  taken  at  Mach  numbers  5  and  8  for  a 
pyramidal  configuration  with  trailing  edge  control  surfaces.  A  sketch  of  the 
test  configuration  is  shown  in  Figure  85*  Calculations  for  the  bottom  surface 
of  this  configuration  were  performed  for  angles  of  attack  of  0  and  10  degrees 
and  control  surface  deflection  angles  of  20  and  40  degrees.  The  average  dis¬ 
tance  of  the  control  surface  hinge  line  from  the  leading  edge  was  used  to 
compute  the  Reynolds  number.  No  attempt  was  made  to  account  for  the  leading 


(43) 


195 


Figure  116.  Comparison  of  Experimental  Data  with  Predicted  Results  for  Blunted  Pyramidal 
Configuration 
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ANGLE  OF  INCIDENCE,  a  , DEGREES 

Figure  116.  Comparison  of  Experimental  Data  with  Predicted  Results  for  Blunted  Pyramidal 
Configuration  (Continued) 
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ANGLE  OF  INCIDENCE,  a  ,  DEGREES 


Figure  116.  Comparison  of  Experimental  Data  with  Predicted  Results  for  Blunted  Pyramidal 
Configuration  (Concluded) 


edge  bluntness  or  the  sweepback.  Figure  116  shows  a  comparison  of  the  analyti¬ 
cal  results  with  experimental  data.  The  comparison  indicates  that  the  calcu¬ 
lated  data  slightly  overpredict  the  control  effectiveness;  however,  in  view 
of  the  fact  that  finite  span  effects  were  not  considered,  generally  good 
agreement  is  shown. 
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Section  VII 


AERODYNAMIC  HEATING  IN  SEPARATED  FLOW 


Flow  separation  brings  about  changes  in  pressure  and  heat  transfer  distribu¬ 
tion.  A  typical  heat  transfer  distribution,  including  the  parameters  used  in 
the  succeeding  correlations,  is  shown. in  Figure  117.  Heat  transfer  rates  de¬ 
pend  on  the  pressure  field.  It  is  well  established  that  surfaces  beneath 
separated  regions  experience  reduced  rates,  while  greatly  increased  local 
rates  are  experienced  in  regions  of  reattachment.  For  a  long  enough  axisym- 
metric  body,  total  integrated  heat  flux  tends  to  approximate  the  flux  that 
would  exist  if  there  were  no  separation  (Ref  730* 


An  investigation  of  the  heat  transfer  in  separated  laminar  and  turbulent  flow 
is  described  in  this  section.  The  model  geometries  considered  are:  a  sharp 
leading  edge  flat  plate  with  a  deflected  trailing  edge,  flap  in  compression, 
and  a  flat  plate  with  an  incident  shock  generated  by  an  external  wedge.  Data 
(obtained  from  various  test  facilities)  along  the  plate  centerline  have  been 
correlated  into  expressions  describing  the  magnitude  and  distribution  of  heat 
transfer  in  the  separated  region. 

DATA  PRESENTATION 

The  local  experimental  heat  transfer  coefficients  were  non-dimensionalized  by 
local  theoretical  heat  transfer  values  for  a  flat  plate  at  zero  incidence, 
evaluated  by  means  of  the  reference  temperature  method. 

Using  the  reference  temperature  relation  (Reference  71 ): 


£-  =  0.5  (J£  +  1)  +  0.044r  M„ 

xoo  T 

oo 

l/2 

where  r  =  Pr  '  =0.85  for  laminar  flow 

l/R 

and  r  =  Pr  '  =0.90  for  turbulent  flow 


The  theoretical  heat  transfer  coefficients  were  obtained  from  (Reference  71) 

1 /p  1/3 

=  Nu*  =  O.332  Re*  '  Pr*  J  for  laminar  flow 
K*  x  x 


and  ~  =  Nu*  =  0,029  Re*  4/^ 
K*  x  x 


for  turbulent  flow. 
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Figure  117.  Definition  of  Interaction  Parameters  for  Shock-Inducod  Separation 


The  downstream  heat  transfer  interaction  length  was  non-dimensionalized  by 
the  attached  boundary  layer  thickness  at  Xq,  defined  by  the  following  simple 
relations  for  laminar  and  turbulent  flow,  respectively: 


-  J. 

5.2  X  Re  2 
o  x 

o 


Laminar 


0.38  X  Re  ~-r-  Turbulent 
o  x  ✓ 
o 


The  evaluated  experimental  data  were  obtained  from  References  7^-80  and 
the  ‘aerodynamic  heating  rates  on  the  test  plates  were  obtained  utilizing 
the  thin  wall  transient  temperature  technique  with  correction  for  heat  con¬ 
duction. 


HEAT  TRANSFER  CORRELATIONS 

It  was  inferred  frcm  experimental  data  that  the  peak  heat  transfer  h  .  and 

the  downstream  heat  transfer  interaction  length  can  be  expressed  by  the 
following  relationships: 


=  f 


Re  ,  (p  , 

.  o  o  , 


=  f 


n  ,.c  ,P  xd] 


o 


P  '  \ 
o  J 


where :  n  =  —  for  laminar  flow 


n  =  *r*  for  turbulent  flow 


Groups  of  dimensionless  parameters  were  plotted  on  log-log  paper  and  the  data 
correlations  were  obtained  by  fairing  a  straight  line  representing  an  average 
through  the  data  distribution.  An  equation  of  the  form  y  -  A  x21,  where  A  is 
the  ordinate- intercept  and  m  is  the  slope  of  the  line,  was  obtained  for  various 
groups  of  parameters,  subsequently  arriving  at  a  correlation  function. 

LAMINAR  SEPARATED  FLOW 


It  was  established  that  the  heat  transfer  rate  decreases  over  the  separated 
region  and  reaches  a  maximum  value  in  the  reattachment  region.  The  magnitude 
of  the  peak  heat  transfer  rate,  location  of  peak  heating,  and  heat  rate  in  tie 
separated  region  are  discussed  in  the  following  sections. 
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Peak  Heating  Rate 


An  attempt  was  made  to  obtain  some  semi-empirical  relations  describing  the 
magnitude  of  the  peak  heat  transfer.  An  evaluation  of  experimental  heat 
transfer  results  leads  -to  the  following  parametric  relation  for  peak  heat 
transfer  (Ref.  8l). 


J*  =  f  (Re, 

fp 

Figures  118  and  119  present  the  peak  heat  transfer  parameter  for  laminar 
separatee  corner  flow  as  a  function  of  pressure  rise  from  the  plateau  value 
to  the  oblique  shock  value  and  from  the  undisturbed  value  to  the  oblique 
shock  value,  respectively,  The  figures  illustrate  that  peak  heating  magni¬ 
tudes  increase  with  decreasing  Mach  number,  and  are  proportional  to  pressure 
ratio.  The  variation  of  peak  heat  transfer  is  easily  obtained  from  Figures 
118  and  119.  In  order  to  determine  the  heat  transfer  dependence  on  Mach  num¬ 
ber,  a  grouped  parameter  (heat  transfer,  Reynolds  number  and  pressure)  is 
shown  plotted  in  Figures  120  and  121  versus  Mach  number.  A  straight  line 
representing  an  average  was  faired  through  the  data,  resulting  in  two  ex¬ 
pressions  for  peak  heating  correlation  functions: 


r&  =  5.8  x  105  I 

[fs 

a0-89  r  **  M  -2.4 

"  1  Rex  Ma> 

(44) 

nfp 

V  P 

/  0 

=  5.0  x  106  | 

h„ 
fp 

(V 

Vpo- 

r  «ex  ^  MJ3-70 

'  O 

(45) 

Downstream  Heat  Transfer  Interaction  Length. 


The  location  of  the  peak  heat  transfer  relative-  to  the  ref&rence  line  is  de¬ 
scribed  by  the  length  H, ^  which  is  illustrated  in  Figure  117.  Figure  122  pre¬ 
sents  the  downstream  heat  transfer  interaction  length  for  M  =  6  and  M  =  8 
conditions  as  a  function  of  pressure  rise.  The  correlation  of  downstream 
heat  transfer  interaction  length  for  laminar  corner  flow  including  the  Mach 
number  dependence  is  presented  in  Figure  123  and  is  given  as  follows: 


(46) 


Heat  Transfer  in  the  Separated  Region 

The  ratio  of  heat  transfer  (hg/hfp)  in  the  separated  region  covering  a  range 

of  Mach  numbers  from  6  to  10  is  shown  in  Figure  124  as  a  function  of  pressure 
rise  from  the  plateau  to  the  inviscid  value.  The  rate  of  heat  transfer  from 
a  separated  laminar  mixing  layer  is  less  than  that  from  an  attached  laminar 
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Figure  118.  Peak  Heat  Transfer  for  Separated  Lominar  Comer  Flow  vs  -s - 1 
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20 6 


Figure  120.  Peck  Heat  Transfer  for  Laminar  Separated  Corner  Flow  vs  Mach  Number 
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Figure  121,  Grouped  Peak  Heat  Transfer  Parameter  for  Laminar  Separated  Corner  Flow 
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Figure  T 23 .  Grouped  Downstream  Interaction  Length  Parameter  for  Laminar 
Separated  Comer  Flow 
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boundary  layer  at  corresponding  values  of  Mach  number,  Reynolds  number,  and 
wall-to- stream  temperature  ratio.  The  ratio  of  heat  flux  in  the  separated 
flow  to  that  in  the  attached  flow  is  a  function  only  of  Prandtl  number  and 
has  been  found  to  be  equal  to  a  value  of  O.56  for  Pr  =  0.72  (Ref-.  82).  An 
average  value  of  heat  transfer  ratio  equal  to  0.425  was  obtained  in  Figure  124. 

TURBULENT  SEPARATED  FLOW 


Peak  Heating  Rate 

Figures  125  and  126  present  data  correlations 
a  function  of  a  characteristic  pressure  rise. 


pressure  ratio 
heat  transfer: 


of  Figure  125  leads  to  the 


for  peak  heating  magnitude  as 
h  k_  V5 

The  plot  of  Re  versus 

il  A  X 

fp  o 

following  expression  for  peak 


P„  0.62 


32.0^ 

0 


N 


(47) 


This  expression  does  not  show  an  explicit  Mach  number  dependence;  however,  a 
Mach  number  influence  is  contained  .Ln  the  P^/Pq  ratio.  A  correlation  of 

h .  P„  -  Pp 

—  as  a  function  of  pressure  parameter  — -p -  shown  in  Figure  126  results 

fp  o 

in  a  relation  for  peak  heat  transfer 

h  .  P_  P_  0.18 

^  =1.18  (-3-§-E)  M„  (48) 

fp  o 


This  expression  does  not  contain  a  Reynolds 
number  effect  is  included  implicitly  in 


number  term;  however,  the  Reynolds 
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Downstream  Heat-Transfer  Interaction  Length 


The  downstream  location  of  the  peak  he|t  transfer  rate  relative  to  the  i;efe 
enee  line  is  determined  by  the  length  jU.  A  parametric  correlation  of  i 2  i 
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shown  in  Figure  127  and  it  is  expressed  by 
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Figure  126,  Peak  Heat  Transfer  for  Turbulent  Fiat  Plate  Flow  with  Extensive  Separation 
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Figure  12/ •  Downstream  Heat  Transfer  interaction  Length  for  Turbulent  Separated  Flow 
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P3/P*  #  OVERALL  PRESSURE  RISE 


Figure  128.  Heat  Transfer  Coefficient  in  the  Separated  Region  for  Turbulent 
Boundary  Layer  -  Shock  Wave  Interaction 
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Heat  Transfer  in  the  Separated  Region 


The  heat  transfer  coefficient  in  the  turbulent  separated  region  would  be  ex¬ 
pected  to  be  higher  than  for  the  laminar  separated  region.  Figure  128  illus¬ 
trates  this  for  data  covering  a  Mach  number  range  from  2.4  to  ?.02.  The 
average  value  of  heat  transfer  for  a  separated  turbulent  boundary  layer  as 
compared  to  an  attached  turbulent  boundary  layer  is  indicated  in  Figure  128 
h 

by  the  ratio  =  0.84. 
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Section  VIII 


CRITICAL  EVALUATION 


CHARACTERISTICS  OF  EXPERIMENTAL  DATA 
Test  Flow  Properties 


The  experimental  data  from  the  several  sources  evaluated  in  this  program  does 
not  completely  cover  the  wide  range  of  possible  local  flow  conditions.  In 
particular,  higher  local  Mach  number  data  for  the  turbulent  cases  would  have 
been  useful. 

For  the  laminar  interaction,  Mach  numbers  varied  generally  from  3*1  through 
8.4,  with  no  sizable  gaps.  There  was  some  data  of  limited  usefulness  around 
Mach  l4.  Local  Reynolds  numbers  ranged  from  1.2  x  105  to  1 .4  x  lO^.  This 
range  is  considered  to  be  insufficient  for  accurately  determining  the  effect 
of  Reynolds  number  variation. 

For  turbulent  interactions,  Mach  numbers  varied  generally  from  1.8  through 
3.8.  There  was  some  data  of  limited  usefulness  around  Mach  6.2.  Local 
Reynolds  numbers  ranged  from  7  x  105  to  6  x  10°.  Neither  the  Mach  number 
variation  nor  the  Reynolds  number  range  of  the  turbulent  data  is  considered 
to  be  wide  enough. 

It  would  oe  desirable  to  have  laminar  separated  flow  data  for  a  test  section 
Mach  number  up  to  between  20  and  25  with  local  Mach  number  reduction  to  be 
achieved  through  pitching  of  the  model.  A  wider  variation  of  Reynolds  number 
may  be  achieved  through  reservoir  pressure  variation  and  variation  of  model 
length . 

Model  Geometry 

In  order  to  simulate  two-dimensional  flow,  high  aspect  ratio  models  are  re¬ 
quired.  The  visual  evidence  of  the  surface  flow  (Section  5-2)  indicates 
that  flow  over  a  finite  span  body  has  three -dimensional  features  and  is 
substantially  influenced  by  the  body  tips .  Thus  the  assumption  of  two- 
dimensionality  seems  questionable  for  all  but  very  high  aspect  ratio  models, 
as  shown  in  Figures  6l  to  69. 

In  some  instances  the  streamwise  spread  of  separated  flow  was  of  the  order  of 
the  model  size  and  no  fully  developed  flow  could  be  obtained .  In  order  to 
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properly  size  the  model,  estimates  of  the  extent  of  separated  flow  should  he 
helpful.  The  length  of  the  forward  part  of  the  model  controls  the  boundary 
layer  thickness,  which  is  a  prominent  interaction  parameter. 

The  addition  of  end  plates  fails  to  eliminate  or  reduce  the  three-dimensional 
effects  and  in  fact  introduces  new  complexities  into  the  flow  picture.  In 
general  the  effects  of  end  plates  on  the  flow  are  not  well  understood  and 
results  are  inconclusive. 

All  three  major  parameters  of  the  flap  geometry  (deflection,  chord  dimension, 
and  span  dimension)  and  flow  variation  should  be  given  detailed  attention  in 
any  future  test  series.  I.t  is  obvious  that  incipient  separation  can  be  pin¬ 
pointed  only  if  flap  deflection  is  varied  systematically  in  small  increments. 
Downstream  Interaction  length  was  found  to  be  sensitive  to  flap  chord  length. 
Flap  effectiveness  can  be  substantially  lowered  by  reducing  flap  aspect 
ratio. 

Flow  Visualization 


Flow  separation  is  a  complex  flow  phenomenon,  and  for  a  thorough  investigation 
some  visual  study  is  essential.  Schlieren  photography  and  oil  flow  surface 
patterns  proved  very  useful  and  should  not  be  overlooked .  There  is  a  need 
for  other  methods  to  indicate  special  flow  details  which  are  difficult  to 
obtain  otherwise  (e.g.  smoke  method). 

Instrumentation 


A  quantitative  picture  of  the  separated  region  is  obtained  by  measurements  of 
pressure  and  temperature.  Two  aspects  concerning  instrumentation  are  of 
paramount  Importance:  instrumentation  arrangement  and  density,  and  measure¬ 
ment  accuracy.  While  availability  of  instrumentation  may  vary  with  the  test 
facility  used,  judicious  arrangement  and  density  of  instrumentation  are  the 
responsibility  of  the  researcher.  Consideration  should  be  given  to  the 
spatial  extent  of  the  phenomena  under  study  and  especially  to  certain  critical 
areas  where  more  flow  details  may  be  desirable. 

ACCURACY"  OF  CORRELATIONS 


Limitations  of  the  data  used  in  constructing  the  correlations  contained  in 
this  report  have  been  mentioned.  The  measurements  have  been  taken  from 
several  sources.  Each  facility  has  its  own  characteristic  peculiarities 
and  generates  its  own  systematic  (as  differentiated  from  random)  errors.  The 
combination  of  data  from  them  all  will  therefore  show  more  scatter  about  the 
best  correlating  line  that  can  be  drawn  through  the  results,  than  is  exhibited 
by  the  data  from  any  single  test  series.  However,  the  resulting  correlation 
is  more  reliable  from  the  inclusion  of  a  large  assortment  of  experimental 
data. 
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Some  remarks  concerning  the  characteristics  of  transitional  cases  will  follow 
shortly.  It  is  stated  here  that  it  is  not  always  possible  to  predict  the 
state  of  the  separating  boundary  layer  and  a  wrong  judgement  on  this  will 
always  lead  to  failure  of  the  correlation  expression. 

Some  of  the  correlations  were  more  difficult  to  construct  than  others. 
Particular  difficulty  was  encountered  in  correlating  the  turbulent  flow  inter¬ 
action  distances.  Therefore  an  idea  of  the  maximum  likely  error  associated 
with  the  use  of  the  formulas  contained  herein  can  be  obtained  by  noting  the 
following.  A  survey  of  the  data  used  in  deriving  the  di/60) turbulent  corre¬ 
lation  showed  that  75  percent  of  the  data  points  had  a  deviation  of  less  than 
+  33%  from  the  correlation  line.  Other  quantities  that  have  been  correlated 
in  this  report  exhibit  much  less  scatter. 

REAL  VEHICLE  CONFIGURATIONS 


It  seems  certain  that  the  data  upon  which  the  results  of  this  report  have 
been  based  do  not  cover  all  flight  conditions  of  interest.  The  working  charts 
which  are  given  here  have  been  constructed  with  this  fact  in  mind,  and  indi¬ 
cate  the  range  of  test  conditions.  It  is  possible  to  select  flight  condi¬ 
tions  that  lead  to  unacceptable  answers,  such  as  d]_  >  Xht.  or  d^  <  1  fi,  and  the 
design  engineer  must  be  prepared  to  exercise  some  discretion  in  using  these 
correlations  (e.g.  Sample  Calculation  Example  2,  Appendix  II). 

Practical  reentry  vehicles  will  exhibit  some  geometrical  features  that  were 
not  reflected  in  the  test  models.  For  example,  surfaces  are  likely  to  be 
highly  swept  as  well  as  blunted,  and  sections  normal  to  the  stream  direction 
may  show  a  steady  taper  in  thickness  from  the  center  to  the  edges.  Such  a 
shaping  will  tend  to  introduce  a  decisively  three-dimensional  character  to  the 
flow.  Figure  84  shows  such  a  configuration  at  an  expansion  angle  of  attack. 

It  is  seen  that  the  flow  is  deflected  in  toward  the  centerline,  away  from  the 
higher  pressure  edges.  Also,  streamlines  are  deflected  about,  either  side  of 
the  rather  low  aspect  ratio  flaps. 

High  speed  vehicles  are  generally  analyzed  on  a  piecewise  basis,  and  this 
viewpoint  should  be  no  less  acceptable  in  interaction  problems.  Conditions 
upstream  of  the  interaction  region  should  be  calculated  as  well  as  possible 
using  standard  methods.  The  correlations  in  this  report  may  then  be 
applied.  Exact  streamline  direction  is  less  important  than  static  pressure 
and  Mach  number.  There  can  be  substantial  reductions  in  flap  effectiveness 
due  to  finite  span,  however. 

SPECIAL  PROBLEMS 


Transition  within  the  Interaction  Region 

In  some  problems,  calculated  transition  distance  will  be  such  as  to  Indicate 
that  the  boundary  layer  at  the  beginning  of  the  interaction  may  be  either 
laminar  or  turbulent.  Transition  nicy  occur  within  the  separated  region.  In 
such  cases  local  pressure  peaks  may  appear.  The  integrated  pressure  is  what 
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is  of  interest,  although  it  is  not  possible  to  say  beforehand  whether  the 
laminar  or  the  turbulent  assumption  is  more,  conservative  (that  is,  leading 
to  a  smaller  integrated  force  increment).  The  turbulent  assumption  leads  to 
higher  local  peak  rates  of  heat  transfer. 

Real  Gas  Effects 


An  experimental  facility  which  produces  the  desired  Mach  number  and  Reynolds 
numoer  for  a  test  model  will  generally  fail  to  produce  the  ratio  of  real  gas 
relaxation  distance  to  characteristic  body  distance  encountered  by  the  real 
vehicle.  Establishing  of  scaling  laws  accounting  for  these  differences  is 
beyond  the  scope  of  this  program.  The  local  outer  flow  properties  should 
be  obtained  from  ideal  or  real  gas  relations  according  to  Section  III. 


Section  IX 


CONCLUSIONS 


Semi-empirical  correlations  describing  pressure  distribution  and  heat  transfer 
in  laminar  and  turbulent  comer  flow  have  been  developed.  These  correlations 
are  based  on  presently  available  experimental  results  covering  a  limited  range 
of  test  conditions.  They  ace  used  for  approximating  the  pressure  distribution 
in  the  separated  region  and  ror  calculating  the  effect  of  flow  separation  on 
aerodynamic  coefficients.  The  developed  prediction  methods  are  applicable  to 
regions  were  the  assumption  of  two-dimensional  flow  is  valid.  In  order  to 
apply  these  methods  to  complex  three-dimensional  configurations,  the  degree 
of  deviation  from  two-dimensional  flow  should  be  established  and  the  proper 
method  of  determining  the  local  flow  properties  selected.  At  some  sector  of 
the  flight  spectrum  the  vehicle  will  encounter  real  gas  effects  which  were 
absent  in  the  test  flows.  These  effects,  which  undoubtedly  will  show  up  in 
the  separated  region,  were  not  included  in  the  available  experimental  data. 

It  is  suggested  however  to  use  real  gas  flow  properties  in  regimes  where  the 
real  gas  effects  become  significant. 

The  primary  emphasis  during  the  development  of  these  correlating  expressions 
has  been  placed  on  providing  the  aerodynamic  designer  with  a  useful  design 
tool  for  predicting  the  aerodynamic  control  characteristics  of  aircraft  and 
reentry  vehicles  operating  at  hypersonic  speeds .  This  method  permits  the 
designer  to  determine  the  presence,  magnitude,  and  extent  of  flow  separation 
and  the  pressure  distributions  induced  by  deflected  control  surfaces.  A 
pressure  distribution  model,  based  on  experimental  observation,  is  employed 
for  determining  force  and  moment  coefficients.  Sample  calculations  are 
included  for  instructive  purposes. 

There  is  some  error  involved  in  the  experimental  measurements  which  were 
evaluated  and  additional  error  is  introduced  due  to  correlating  procedure. 
Extrapolation  of  the  prediction  expressions  outside  the  range  of  test  con¬ 
ditions  is  feasible  but  should  be  accomplished  with  caution. 


Appendix  I 


INCREMENTAL  FORCE  AND  MOMENT  EQUATIONS 


The  following  is  a  derivation  of  the  equations  for  the  two-dimensional  force 
and  moment  increments  produced  by  separated  compression  corner  flow.  The 
idealized  pressure  distribution,  coordinate  axes  system,  and  force  and  moment 
convention  used  in  the  derivation  are  defined  in  Figure  129. 
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2.  Pitching  moment  increment  about  !IL  due  to  normal  force  increment: 
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3.  Center  of  pressure  of  normal  force  increment: 
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5.  Pitching  moment  increment  about  HL  due  to  axial  force  increment: 
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6.  Center  of  pressure  of  axial  force  increment: 
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7.  Total  increment  of  pitching  moment  about  MRP: 
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8.  Windward  component  of  elevon  hinge  moment: 
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Appendix  II 
SAMPLE  CALCULATIONS 


The  solutions  to  two  sample  problems  are  worked  out  in  detail  to  illustrate 
the  procedure  for  calculating  control  effectiveness  which  is  developed  in 
Section  6.  A  simple  hinged,  flat-plate  configuration  is  analyzed  for  separate 
flight  conditions.  Example  1  is  a  straightforward  problem  for  which  the  me¬ 
thod  is  directly  applicable.  It  is  shown  in  Example  2  that  for  certain  flight 
conditions  the  separated  region  extends  forward  to  the  leading  edge  of  the 
configuration.  For  this  case  a  modification  of  the  method  is  suggested  which 
provides  a  solution. 

Example  1. 

Given: 

0  2  4  (ft)  6  8  10 


I  (  >  t  i  i 


Configuration:  Two-dimensional  flat  plate  with  negligible 
thickness.'  Span  is  b  ft 


Moo 

10 

’Altitude 

= 

150,000  ft 

^co/  ft 

= 

1.06  x  105 

Poo 

= 

2.84  psf 

T* 

= 

479' “R 

a 

ss 

15  deg. 

T 

wall 

= 

1500 °R 
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Problem:  Calculate  the  incremental  normal  force  and  pitching 

moment  coefficients  produced  on  the  windward  surface  by  a 
positive  flap  deflection  of  ten  degrees.  (g^  =  +10°). 

Assume  negligible  viscous  interaction. 


Step  1 .  Determine  the  local  flow  conditions  . 


P 

a 

Poo 


=  13.5 


(Figure  (96c) 


^Cp  Ja 
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£  Moo 


=  0.179 


Toe 


Re. 


a 


Re. 


3.0 

(Figure  97c) 

1440  °R 

5.5 

(Figure  98c  ) 

1.94 

(Figure  99c) 

00 


—  =  2.06  x  nr 

ft 

Step  2.  Verify  the  existence  or  non-existence  of  separated  flow. 

(0  )2  =  0.11  (Oblique  shock  compression) 

inviscid 

Rea  =  8  (2.06  X  105)  =  1.65  X  106 

HL 

For  simplicity,  assume  Rea  transition  =  2  X  10^  (For  a  discussion  of 

transition  Reynolds  number,  see  Section  4.1) 


C  =  0.020  (Figure  104) 

P« 

incipient 


Therefore,  since 


>  C 

incipient 


separation  does  exist. 
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Step  3 .  Determine  where  the  flow  separates . 

For  the  1st  iteration^  assume  X  1  =  7  ft 

o 

Re  =7  (2.06  x  105)  =  1.44  x  106 


(c  )  =  0.0185 

p«  1 1 


(Figure  106) 
(Figure  108) 


T 

wall  _  , 

m"  •** 


X 

Re* 

a 

Rea 

= 

0.28 

(Figures  100  &  10l) 
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ft 
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n lie 
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0.0215 

(Figure  102) 
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O.057  ft 

^i\ 

= 

1.14  ft 

i 

(O 

1 

= 

8.14  ft>X’ T 

iLL 

Assume  X^  =  6.5  ft  for  2nd  iteration. 

Rea  =  1.34  x  106 
X’ 


0.019 

21.5 

0.055  ft 
1.18  ft 
7.68  ft  <  X* 

rLb 
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X1  +  d 
o 


X‘ 

o 


7.5 


Step  4.  Calculate  the  windward  pressure  distribution 
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(Figure  110 ) 


(Figure  112) 
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Equivalent  flow 


deflection  angle  for(c  )  = 

pap 

2.6  deg  (obtained  from  oblique 
0.35 

0.42  ft 


0.019: 

shock  relations) 

(Figure  114) 


(C  )_  =  0.11  (oblique  shock  compression) 

p  inviscid. 

(C  )p  =  0.114  (compressing  isentropically) 
p  inviscid 

Since  <  15°  >  use  the  isentropic  value .  (See  Section  IV) 

/M  \2  P 

(C  >2  =  Iff)  F  +  (Ct>  } 

pco  2  pa  2  PCD  a 

=  0.644 


X'  =  -  (X-5) 


Step  5 •  Evaluate  the  incremental  forces  and  moments . 
1.  Normal  force  increment: 
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Ac, 


i,sL 


=  [0.257  -  0.179]  [  1.2  -  — ^  +  0.985(2)] 
+  [o.644  -  0.257]  (0.985)  [2  - 


=  0 .796  ft 

2.  Increment  of  pitching  moment  about  HL  due  to  normal  force  increment: 
AC_  SLj 


lnHL  ^ANw 


<<Wp  -  'V* 


rift  4idi ,  di2  V?' 

6  '  2  2  2 


rr  2  2 

d  2] 

(Vs  ■  (  Vp 

(cos2  8f } 

Cf  d2  d2d3 

2  “  0  , g 

3 

"  “TT 

-  [0.257  -  0.179] 


(0.45)2  (Q.45)(1.2)  (1.2)2  (o.985)2(2)2 

6  "  2  2  "  2 


-  [0 . 644  -  O.257]  (0.985)* 
_  (0.42)21 

o 


l£i  .  _(%6ir 
2  6  6 


=  -0.813  ft 

3.  Center  of  pressure  of  normal  force  increment: 


AC  L 
nHL 


m,„  '] 


An 


°pAn 


-  XH1,+ 


w 


Ac 


N 


„  -0.813 

=  ”3  +  ~o?f9& 


w 


=  -4.02  ft 


4.  Axial  force  increment: 


Ac.s 

A 


I  v  /  _ .  (r. 

*  n  _  V  I  III 


=  (.smd. 


+  teJ*-(cJp]  [°* '  r  -  ?\\ 

=  (0.174)  | [0.257  -  0.179)  2  +  [o.644  -  0.257]  [2  -  -~i  -  ~=Jj 


=  0.127  ft 
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5.  Increment  of  pitching  moment  about  HL  due  to  axial  force  increment: 


AC  SL 
mHL  JAA 


-  =-(^8f)j[(cpJp-(cJK]f/  + 
2  d22  d2d3 


^Wp- 


2  &  6  b  j  j 

=  -  (0.174)2  J[o. 257-0. 179]  +  [o.644-0.257] 

L%f  _  (0-61)2  .  (0.61)^0.42)  _  (o^42)2j| 


=  .-0.026  ft 

6.  Center  of  pressure  of  axial  force  increment: 

AC  Ll 

mHL  J. 


AA 


W  zz  7  -  _ 

cp  HL  AC, 


v  -Q  .026 

—  m  0  - 


0.127 


w 


=  0 .204  ft 

7.  Increment  of  pitching  moment  about  MRP: 


VI  A_!ix  .  %lw 

b  b  CP&N 


w 


.  z  =0.796  (-4.02) 
b  0p4A 


w 


2 


-  0.127  (0.0204)  =  -3.23  ft 

8.  Windward  component  of  elevon  hinge  moment: 

2 


HM 

bq 


Jw 


"  •  j(Cpj2"2  '  Z  [^pJs-^pjJ  (d22  +  d2d3  +  d32)] 

=  -  jo. 644  -^|^2 -  \  [0  .644-0.257]  [(0.61)2  +  (o.6i)(o.4a)  +  (o.42)$ 
=  -  1.24  ft2 
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Configuration:  Two-dimensional  flat  plate  with  negligible 
thickness . 


M 

00 

= 

20 

Altitude 

= 

200,000  ft 

Re®/ft 

= 

3.27  X  10 ^ 

Poc 

= 

4.13  X  10“1  psf 

= 

457  °R 

a 

= 

10  deg. 

Twall 

s 

2000  °R 

Problem:  Calculate  the  incremental  force  and  pitching  moment  coef¬ 
ficients  produced  on  the  windward  surface  by  a  positive 
flap  deflection  of  ten  degrees  (Sf  =  +10°).  Assume 

negligible  viscous  interaction. 

Step  1.  Determine  the  local  flow  conditions. 


8°li? 

II 

21 

£*  .  ! 

(Figure  96d) 

:  ) 

>--0-  =  0.071 
l  Moo 

11 

4.3 

(Figure  97d) 

T«  = 

1970  °R 

M  a  = 

9.8 

(Figure  98d) 

1.87 

(Figure  99<i) 

6.11  X  10 


4 
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Step  2.  Verify  the  existence  or  non-existence  of  separated  flow. 

(C  )5  =  O.O87  (Oblique  shock  compression) 

Pa 

Re. 


'a 


=  8  (6.11  x  101*)  =  4.9  x  105 


HL 

Therefore  the  flow  is  laminar.  (See  Figure  24) 

C  =  0.019 

Pa 

incipient 

Therefore,  since  (C  )0  >  C  ,  separation  does  exist, 

^  ^incipient 

Step  3*  Determine  where  the  flow  separates. 

For  1st  iteration,  assume  X'  =  7  ft 

°  4  5 

=  7  (6.11  X  10  )  =  4.28  X  10? 


(Figure  104) 


Rea 

=  7  (6.: 

X' 

°1. 

(c  ) 

*  0.020 

Pa  Pi 

/  j  \ 

ft 

=  l4o 

(Figure  106) 
(Figure  108) 


T 


wall 


La 


~  1 


* 

Rea 

Re^, 

=  0.079 

(Figures  100  and  101 ) 

a 

* 

^a 

=  4.83  x  103 

ft 

8 

=  0.075 

(Figure  103) 

-Jr 

=  0.199  ft 

(di)1 

=  27.9  ft 

(di) 

x  1 

=  34.9  >xHL 
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Assume  X' 
o 


(C  ) 


P«  P2 


4  ft  for  2nd  iteration 
4  (6.11  X  104)  =  2.45  X  105 

0.023 
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0.15 

34.5 

38*5  >XHL 


Assume  X' 
o 


(C  ) 


P«P3 


2  ft  for  3rd  iteration 
2  (6.11  X  104)  =  1.22  X  105 

0.02? 


(8°l 

=  4oo 

8°3 

=  0.106 

(a.). 

'  X'3 

=  42,4 

♦  <v3 

-  kk’li  >  XHL 
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XU1  *  X'  +d 
HL  o 


The  magnitude  of  indicates  that  realistically  the  region  of  separation 
extends  forward  to  the  leading  edge  of  the  plate,  for  this  case,  an  alternate 
approach  can  be  used  which  is  based  on  the  assumption  that  the  plateau  pre- 
sure  region  extends  forward  to  the  leading  edge;  i.e.,  d^  =  and  =  0. 

Step  4,  Alternate  method  for  calculating  pressure  distribution. 

Let  d_  =  X’  =  8  ft 

I  nJLi 


=  C.25 


<  0.25 


2  ~ ■  =  0.475 

dl 


(Figure  112) 


i 


c 


ft 


<*>= 


tan*"'1' 

tan"'1' 

.  -1 
tan 


d2  sinSf 
ai  +  d2  cos8f 

1.9  sin  10° 

8  +  1.9  cos  10s 

0.0335 


2° 


Assuming  the  plateau  pressure  corresponds  to  wedge  pressure  for  a  local 
flow  deflection  angle  of  4>  degrees, 


However, 


(V;p 

=  0.0086  (Obtained 

from,  oblique  shock  relations) 

/ 

M  \2 

P 

^cpjp 

-  u( 

t) 

a 

+  (C  ) 

P®  a 

=  0.0086  ^ 

9.8f 

20  i 

21 

+  0.071 

=  0.114 

'  / 

Let  ifi  «  0 
d 

=  0.25  (Figure  114) 

dl 

d3  =  2  ft 


since  d^  must 

be  less  than 

or 

equal  to  dp,  let 

dj  =  d2  = 

1.9  ft 

it 

CM 

'~h 

o°< 

V — ^ 

0.087  (Oblique  shock  compression) 

/  M  f 

P 

^Cpco^2  = 

(op“^  ) 

a 

+  (C  ) 

N  p®  « 

21  + 

0.071 

0.51 

X*  =  -  (X-5) 

Step  5.  Evaluate  the  incremental  forces  and  moments. 
1.  Normal  force  increment: 


[(C  )  ~(C  )  1  f  d.  - 

|  P®  P  p®  a  j  i  1 


ifi 


’cosSf  Cf 


c  -  ^ 

Of  2 


[o.ll4  -  0.071]  [8  +  cos  10°  (2)] 

+  Jo. 510  -  O.IIU]  (cos  10°j[  2  -  ^  ] 

0.468  ft 
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2. 


Increment  of  pitching  moment  about  HL  due  to  normal  force  increment: 


Ac  SL 
m 
HL 


,An 

J  w 


(c  )  -(c  ) 

Poo  p  Poo  a 


A i2  if'A .  “A  cf2’ 
“5“  "  2  +  2  2 


2  2  2 
cf  d2  d2d3  d3 


-  (Cpj2-(CpJp]<°0s28f> 

=  [o.n4  -  0.071]  j 


z — z~  -  r 


-  [0.510  -  0.114]  (-cos2  10°) 
=  1.215  ft2 

3.  Center  of  pressure  of  normal  force  increment: 


°pAn 

w 


M2. 

2  6^  6^ 


Ac  Ll 

4^ 

'HL  AC. 


N 


w 


„  .  1.215 
=  '  3  +  o^5S 


=  -  0.4  ft 


4.  Axial  force  increment: 
Acas 


c  -fs .  !al! 

f  2  2  jj 


=(sinSf)[[(cpJp-(cpjQj  °f  +  [<s„)2-(cpJp] 

=  (0.174)  j [0 . 114  -  0.07l|  2  +  [0.510  -  0.114 J  [2  -  -  -~2j| 


=  0.022  ft 
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<s. 


5.  Increment  of  pitching  moment  about  HL  due  to  axial  force  increment: 


AC  SL 
m 

HL 


AA  o  (r  C  „2 

-  -  -KM  KV(C».4]  -i  +  KV(cJp] 


2  ?  2 

cf  d2-  d2d3  d3 

2 ~~z - 5"  “ r 


=  -  (0.174)2  ([0.114-0.071]  +  [0.510-0.114] 

"ill2  U;9)2  U;9)2  (1.9)2^ 

2  ST”  "  6  S~ 


=  -  0.005  ft 

6.  Center  of  pressure  of  axial  force  increment: 


Z  =  Z 

cpAA 
w 


HL  Ac 


Ac  LI 

mHL 

AAw 

w 


=  0  - 


-0.005 

0.022 


=  0.23  ft 

7.  Increment  of  pitching  moment  about  MRP: 


Ac  SL  A  n  a  1 

fnMRP  Jw  _  ^CMS|w 


X 


cpAn 


w 


Ac.s  . 

A  w 


CpAA 


w 


A  l./^O  /  A  ).  \  A  AAA  /  A  AA\ 

=  u.hOu  -  \j  *\jdc 1 

2 


8. 


=  -O.192  ft 

Windward  component  of  elevon  hinge  moment: 


HI*  . 

bq 


-  H  )2  ,c£_  . 

■  Pco  *=  2 


t  ^  ld22  +  d2d3  +  d3S] ! 

12 

0.510  &  -  z  (°* 510-0.114)  {(1.9)2  +  (1.9)2  +  (1.9)2 


=  -0.306  ft£ 
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